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Abstract. Starting from hyperbolic dispersion relations, we derive a system of Roy-Steiner equations for 
pion Compton scattering that respects analyticity, unitarity, gauge invariance, and crossing symmetry. 
It thus maintains all symmetries of the underlying quantum field theory. To suppress the dependence of 
observables on high-energy input, we also consider once- and twice-subtracted versions of the equations, 
and identify the subtraction constants with dipole and quadrupole pion polarizabilities. Based on the 
assumption of Mandelstam analyticity, we determine the kinematic range in which the equations are valid. 
As an application, we consider the resolution of the 77 — > mr partial waves by a Muskhelishvili-Omnes 
representation with finite matching point. We find a sum rule for the isospin-two S'-wave, which, together 
with chiral constraints, produces an improved prediction for the charged-pion quadrupole polarizability 



(Q2 - P2) 7 ' 



(15.3 ± 3.7) • 10 4 fm 5 . We investigate the prediction of our dispersion relations for the 



two-photon coupling of the cr-resonance -T CT77 . The twice-subtracted version predicts a correlation between 
this width and the isospin-zero pion polarizabilities, which is largely independent of the high-energy input 
used in the equations. Using this correlation, the chiral perturbation theory results for pion polarizabilities, 
and our new sum rule, we find r all = (1.7 ± 0.4) keV. 



PACS. 11.55.Fv Dispersion relations 
14. 40. Be Light mesons 



11.80.Et Partial-wave analysis - 13.60.Fz Compton scattering 
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1 Introduction 



The reaction 77 — > tttt is of particular interest in the realm 
of non-perturbative QCD, as it probes strong-interaction 
dynamics in the ++ channel, which has the same quan- 
tum numbers as the QCD vacuum. However, a theoretical 
understanding of the dynamics in this channel has long 
proven elusive. The pions produced in the fusion of two 
photons are strongly interacting, such that a description of 
the reaction 7T7t — > tttt is a prerequisite for understanding 
77 — > 7T7T. A significant step forward in this direction was 
made in [1] , where it was shown how combining constraints 
from the analyticity, unitarity, and crossing symmetry of 
relativistic quantum field theory (the Roy equations [2]) 
with the chiral symmetry of QCD (using Chiral Perturba- 
tion Theory, ChPT [3111]) produced detailed information 
on pion-pion amplitudes in this channel. As a result a 
very precise prediction of the pole mass m a = M a — iT a /2 
of the <7 resonance — the lowest-lying resonance in QCD — 
was obtained 



441^ 6 MeV, 



544^ MeV. 



(1) 



This resonance influences the cross sections in 77 — > 
7T7T, which therefore provides an alternative to meson- 
meson scattering reactions for its excitation. Experimen- 
tally, 77 — > tttt is accessible in e + e~ colliders via the re- 
action e + e~ — > e + e~7T7r, where both the incoming elec- 
tron and positron radiate one photon [5l[6l[7ll8ll9lfT0] . How- 
ever, due to its very large width, the a is only manifest 
as a broad bump in the 77 —> 7r°7r° cross section, which 
makes it difficult to extract information on the resonance 
from these data. Therefore, an improved theoretical un- 
derstanding of this important reaction that respects all 
available constraints is strongly called for. 

Moreover, the two-photon coupling of the a influences 
Compton scattering from the proton, via the possibility 
to exchange degrees of freedom corresponding to the res- 
onance between the incoming photon and the target pro- 
ton [TTlfPTlfRTlIT^] . Developing a theory of 77 — > tttt that 
includes the ++ resonance, and delineates its influence 
on cross-section data in this channel, is therefore an im- 
portant step. 

Besides its relation to nucleon polarizabilities, there 
has been particular interest in the two-photon width of the 
a as inferences are then made concerning the nature of this 



2 



M. Hoferichter, D. R. Phillips, C. Schat: Roy-Steiner equations for 77 — s> tttt 



state. Apart from an interpretation as a qq state, possibil- 
ities such as a tetraquark state, a meson-meson molecule, 
or a gluonic resonance have been suggested in the litera- 
ture. The coupling to two photons has been used to dis- 
criminate between different scenarios (see [T51IT5] and ref- 
erences therein). Extraction of the cr's two-photon width 
F ail from data on 77 — > tttt thus becomes an important 
piece of this puzzle. A recent X-matrix approach to the 
extraction of a widths from 77 — > tttt data can be found 
in [T6lfl7] . Alternatively, model-independent 77 — > tttt dis- 
persion relations have been obtained and solved by means 
of Omnes techniques [T5UT§] in descriptions of this reac- 
tion [^[^[^[^[Ml[^[T^[Ml[^[^[ni[29l[3n] . The most 
sophisticated such treatment was in [315], where, motivated 
by the fact that most of the Belle data lie above 1 GeV, 
a Muskhelishvili-Omnes representation was constructed 
that dynamically includes the KK channel, in order to 
obtain a description of 77 — > tttt valid up to 1.3 GeV. 

In this work, we consider a more general approach, 
namely a complete system of Roy-Steiner equations for 
77 — > tttt and the crossed reaction ^tt — > 771 that, in anal- 
ogy to the 7T7r Roy equations, fully respects analyticity, 
unitarity, and crossing symmetry of the scattering ampli- 
tude. We find that — at a similar level of rigor at which 
the Roy equations for tttt scattering hold — our 77 — > tttt 
equations are valid up to 1 GeV. (The domain of validity 
can be extended further under certain additional assump- 
tions.) We compare our equations for the 77 — > tttt par- 
tial waves to existing dispersive descriptions of this pro- 
cess, and find, in particular, that a numerically important 
coupling between S- and D-waves has been previously ne- 
glected. Furthermore, our equations lead to sum rules for 
the isospin-two partial waves, which we use to improve the 
ChPT prediction of the charged-pion quadrupole polariz- 
ability. 

As an application of these results, we study the con- 
straints of our equations on the two-photon coupling of 
the a. The subtraction constants necessary to ensure suf- 
ficiently fast convergence of the dispersion integrals can be 
directly related to the pion polarizabilities, which there- 
fore play a similar role to that of the tttt scattering lengths 
in the case of tttt Roy equations. As the tension between 
various experimental determinations of the dipole polariz- 
ability of the charged pion (based on Primakov measure- 
ments [31] or radiative pion production [321133] 1 and ChPT 
predictions [34 , 35 ( 36 ( 37, 38 , 39 is far from being resolved, 
we provide the two-photon width of the a as a function 
of the pertinent polarizabilities. We are confident that the 
ongoing measurements at COMPASS 10 11. in combina- 
tion with ChPT predictions, will clarify the situation in 
the near future. 

The paper is organized as follows: having specified our 
conventions in Sect. [3J we present the detailed derivation 
of the Roy-Steiner system in Sect. [3] The domain of valid- 
ity of these equations is studied in Sect. SJ We then con- 
centrate on the equations for the 77 — > tttt partial waves, 
whose solution in terms of Omnes functions is discussed 
in Sect. [5j Establishing the connection to the two-photon 
width of the a in Sect. [5] we then discuss the input we 



use and present our numerical results in Sects. [7] and [5] 
Various details of the calculation are relegated to the ap- 
pendices. 

2 Formalism 

2.1 Kinematics 

We first consider the Compton-scattering process 

7(«i, Ai)7r tt (pi) -> 7(92, A 2 )7r h (p 2 ) (2) 

with momenta as indicated, photon helicities Ai, A2, and 
pion isospin indices a, b. For on-shcll particles, the Man- 
delstam variables defined as 

s=(Pl+qi) 2 -, t=(q 1 -q 2 ) 2 , u = (q 1 -p 2 ) 2 , (3) 

are subject to the constraint 

s + t + u = 2Al 2 . (4) 

In the center-of-mass frame (CMS), we have 

t = -2q 2 (l-z s ), z s = cos6 s , q 2 = [ - (5) 

4s 

The S'-matrix for the charged process can be written as 

out (7(92, A 2 )7r ± (p 2 )|7(9i, Ai)7r ± (pi)) in 
= (27r) 4 <5 4 (q 2 +P2 ~ qi - Pi) 
x{d XlX2+i e 2 F c XiX2 ( S ,t)e^-^}, (6) 

with e 2 = Ana and azimuthal angle ip, and similarly for 
the neutral amplitude F™ ^ . Separating the photon po- 
larization vectorfQ, 

F Xl x a («, t) = Xi)el(q 2 , \ 2 )W^(s, t), (7) 

we can use gauge and Lorentz invariance to decompose 
the amplitude as 

W^(s,t) = A(s,t)(^g^ + q 2 ^qiu) + B{s,t)(2tA f ,A l/ 

- (s - ufg^ + 2(s - ?i)(Z\ M (ji„ + A v q 2l Sj, 

(8) 

where = pi^ + p 2ll and we have dropped terms that 
vanish in F Xl \ 2 due to e(qi, Ai) • qi = 0. In the conventions 
of [12] for the polarization vectors, one obtains 

F ++ (s,t) = F__(a,t) = 4(Af 4 - su)B(s,t), (9) 

- - ~A(s,t)+t(t-AM 2 )B(s,t). 

1 Here and below, we suppress isospin indices whenever pos- 
sible. 
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The Born-term contributions are 
A BoTn (s,t) ee A c ^ Bmn (s,t) 

= 2tB c ' Bom (s,t) ee 2tB amn (s,t). (10) 



M 2 - s M*-u 

jBc 



In these conventions, the differential cross section is given 

by 

^ = ^-(\F ++ ( s ,t)\ 2 + \F + _( s ,t)\ 2 ). (11) 

The analytic continuation of F + ±(s, t) to the kinemat- 
ical region where t > AMI describes the crossed-channel 
process 77 — > mr, such that the cross section reads 



da 
di? 



^-a(t)(\Fl + (s,t)\ 2 + \Fl_(s,t)\ 2 ), (12) 



where a(t) = yjl - 4M%/t. The formula is the same, al- 
beit with an additional factor of 1/2 on the right-hand 
side, for the case of neutral pions. It should be noted, 
though, that the kinematics for the crossed reaction 



7(91, Ai)7(-g 2 , Aa) -> TT a (-p 1 )TT b (p 2 ) 



(13) 



in terms of the Mandelstam variables (jSJ) lead to different 
center-of-mass momenta for initial and final states 



2 t 
q 2 = P 



and to the CMS scattering angle 
v 



z t = cos f t 



Aptqt 



V = S-U, pt = I Pi 



1i 



(14) 



|q*|. 

(15) 



2.2 Partial-wave expansion and pion polarizabilities 

The partial-wave expansion of the amplitudes for pion 
Compton scattering reads [4"5] 

00 

F+±(s,t) = £(2J+1)/ Jl ±(a)df j±1 (2!,) ) (16) 
7=1 

with the Wigner d-functions0 



j(j+xy 



and the inversion 
1 



fj,±{ s ) = \f dz s dl ±1 (z s )F +± {s,t) 



i=-2q 2 (l-z 3 ) 



(17) 



(18) 



The expansion (|16[) can be mapped onto the multipole 
expansion [H] vi£0 



/J,±(*)=± 



2^i 



a(2J+ 1) 



(-Bj(w) ± Mj(w)), (19) 



2 For convenience, we write m ; (cos 8) instead of d J m m , (8). 

3 The covariant amplitudes in [44] are related to ours by 



where ui denotes the energy of the photon. Defining the 
pion polarizabilities as the leading terms of the Born-term- 
subtracted multipoles Ej{uj), Mj(uS) in an expansion in uj 

EH, 



otj 



2J[(2J-1)!!] 2 Ej( 



UJ 



J+l 



,2J 



2J[(2J- l)!!] 2 Mj( 



J + l 



UJ 



2,1 



(20) 



j=0 



we can read off aj and j3j from an expansion of the Born- 
term-subtracted amplitudes F + ±(s, t) in t at fixed s = Ml 

2a - t 
— F ++ (s = Mlt) ^ ai +f3 1 + -{a 2 + ft) + 0(t 2 ), 

^*+-(« = Mlt) = ai -0 1 + ±( a2 - ft) + 0(t 2 ). 

(21) 



2.3 Relation to 77 — > tttt 

To establish connection to the notation of the crossed pro- 
cess [3"UII38| . which we will refer to as the t-channel reac- 
tion, we briefly discuss 

7(gi,A 1 ) 7 (g 2 ,A 2 )^^(p 1 )7r b ( ? , 2 ) (22) 

in terms of the Mandelstam variables 

S = (qi+q 2 ) 2 , t = {q 1 -p 1 ) 2 1 u = { qi - p 2 ) 2 , (23) 

and the amplitudes 

out(7r(pi)7r(P2)|7(gi> Xih{q 2 , A 2 )) in (24) 

= ie 2 (2n) i S i (q 2 +p 2 - q x - Pl )H Xl x , 2 (s, t)e^-^ . 

They are related to the s-channel amplitudes by 

#++(*,*) = -F + _(t,s), H + _(s,t) = -F ++ (t,s), (25) 

and their polarizability expansion therefore reads 

2a - s 
—H ++ (s, t = Ml) = ai - ft + — (a a - /3 2 ) + 0(S 2 ), 
M n s 12 

—2a - s 

—H+-(Z,t = Ml) = ai + ft + — (a a + A) + 0(§ 2 ). 

Mas 12 

(26) 

Furthermore, the partial- wave amplitudes hj± follow from 
F + ± via 

F ++ (s,t) = -^(2J+l)^_(<)rf 2 J (z t ), 
J 

F + _(s,t) = -J2( 2J + l)hj, + (t)4oo(*t), (27) 
j 

where due to Bose symmetry the sum extends over even 
values of J only, and 



di (z) = Pj(z), di (z) 



Agk = -e J A, B G r = 16e A B. 



2P! J _ 1 (z)-J(J-l)Pj(z) 

y/(J-l)J(J+l)(J + 2) ' 

(28) 
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In our conventions, the transition between isospin and par- 
ticle basis is achieved by 



(29) 



3 Roy-Steiner equations 

First we review the basic steps in the construction of Roy 
equations in tttt scattering. They are [2J: 

1. Write down a twice-subtracted dispersion relation at 
fixed Mandelstam t, whose subtraction "constants" are 
actually functions which depend on the value of t cho- 
sen. 

2. Use crossing symmetry to determine the subtraction 
functions, such that the remaining free parameters are 
the tttt scattering lengths. 

3. Expand the imaginary part of the amplitude that ap- 
pears under the dispersion integrals in partial waves, 
and perform a partial-wave projection of the resulting 
equation. 

In this way, one arrives at a system of integral equations 
for the 7T7r amplitudes tf(s) 



n OO 

2 oc p 

#)+EE / ds ' K iv 



(s,s')Im#V)> ( 30 ) 



that relates a partial wave of a given angular momentum I 
and isospin / to all other partial waves via analytically cal- 
culable kinematic kernel functions K(p (s,s'). A detailed 
discussion of how to numerically solve this system can be 
found in |45j . 

The construction of Roy equations for reactions with 
non-identical particles is hampered by the fact that cross- 
ing symmetry intertwines different physical processes, 
such that the second step based solely on fixed-i dispersion 
relations fails. For this reason, hyperbolic dispersion rela- 
tions [46j were used in [471148] to determine the subtraction 
functions and derive Roy equations for ttK scattering, as 
within the hyperbolic approach the dispersion relations 
automatically involve both ttK — > ttK and 7T7t — » KK 
physical regions. The resulting system in such a case is 
therefore referred to as Roy-Steiner equations. 

In this work, we will go a step further and solely con- 
sider hyperbolic dispersion relations. This is particularly 
convenient because s ■<-> u crossing symmetry is manifest, 
and both 771 — > jtt and 77 — ¥ tttt amplitudes contribute, 
such that all constraints by crossing symmetry are auto- 
matically fulfilled. We now turn to the setup of this system 
of Roy-Steiner equations. 



3.1 Hyperbolic dispersion relations 

We start by writing down unsubtracted hyperbolic dis- 
persion relation for the amplitudes A and B, which can 



be constructed following [15] . The advantage of using dis- 
persion relations for A and B instead of F + ± is that all 
constraints by gauge invariance that lead to the decom- 
position (0) are automatically built in. In particular, the 
equations for F ++ and F_| do not decouple, as gauge in- 
variance dictates that the same invariant function B con- 
tributes to both amplitudes. The dispersion relations for 
A and B read 

00 

A(s,t)=A B (s,t)+ 1 - j dt ' l ^^l 

AMI 



00 

i / ds'lmA(s',t')(^- 

Ml 



s s — u s — a 



, ^ RH , , 1 f , ,lmB(t',z') 
B(s,t)=B B (s,t) + - / At' 



TT J t'-t 
4M? 



(31) 



OO 

i J ds'lm ^73 



s s — u s — a 



with the Born terms 



A B (s,t) =A Bmn (s,t) 



1 



Ml -a 



B B (s,t)=B Bo ™(s,t) 
(32) 

only contributing to the charged-pion process (^4 Born and 
S Born were already defined in (flO|) '). The primed set of 
Mandelstam variables is constrained to lie on the hyper- 
bola 

(s — a)(u — a) = (s — a)(u — a) = b, (33) 

where the hyperbola parameter a can be freely chosen. In 
particular, it can be used to optimize the range of va- 
lidity of the resulting system of Roy-Steiner equations 
(cf. Sect. [4]). The above integrals are understood such that 
the integrands shall be expressed in terms of the integra- 
tion variable and the external kinematics by virtue of Q33p 
and 

s'+t' + v! = 2M*. (34) 

The second integral in pip is reminiscent of fixed-i disper- 
sion relations, but in that case JmA(s',t') —> lmA(s',t) 
and the last term is removed. Thus, the key difference here 
is that t' depends not only on t, but on s and s' as well. 
However, it is possible to recover the limit of fixed-t by 
sending a to infinity. This can be shown explicitly based 
on the relation between z s and z' s given in (|A.6j) . 



3.2 Sum rules and subtracted dispersion relations 

The most economical way to obtain a subtracted version of 
hyperbolic dispersion relations is to derive sum rules from 
the original dispersion relation and then subtract them 
from it. We will choose the subtraction points such that 
the subtraction constants coincide with the pion polari- 
zabilities. 
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With knowledge of the dipole pion polarizabilities ai± functions for the unsubtracted case read 



Pi we can implement one subtraction, while for a second 1 

subtraction the quadrupole polarizabilities a2±/?2 are also ++ f sq 2 2J' + 1 f j 

needed. For example, choosing s = M 2 and taking the Kjj'{ s ' s ) = s /q/2 2 J dZs ^ Zs )^ii( Zs ) 

limit t — ► 0, we can compare (f5T|) with (0) and ([2~Tj) in -1 

order to obtain^ ,j> , , 



x 



M-K 4M 2 /" ,ImB t' z L 
^L( ai+ ^ 1 ) = ^L / di' 1 ' tJl ° (35) 



1 + z' I s' — s s' — u s' — a 
1 



4Ml 

00 

4M[ 



G ^ (S '* ,} = ^ 4A/ 2 ) 2 / dZs(1 + ^iiW 



^ / d S 'ImS( S ',O| (- 



M 2 s'-cj' rf 20 



1 

,2 



which, together with similar sum rules, may be used to — \-r„ s '\ _ q 2 J + 1 /" [\ — z W J f z ) 
write down subtracted versions of (EU). JJ '^ ' ' 4s'q' 2 2 J sl lJ sJ 

-1 

x 11V s ' (t' - 1 L 

1 + [ s — s s — U s' — a J 



3.3 s-channel projection 



1 

.2 o 7' 



, q 2 J + 1 f j 

The Roy-Steiner system is obtained by expanding the ^jj'( s ^ s ) — ~ 2 / dZs ^ ~ z s)di,-i( z s) 

integrands into partial waves and subsequently project- 
ing each equation onto s- and i-channel partial waves. To 
achieve this we use (1T61) and (l27ll as well as their inver- 



-1 



di'M r 1 1 1 



sions. It is useful to note that 1 — z ' s I s — s s — u s — a 

1 

4(M* - au) = 8 S q 2 (l + *.) = -t(t - 4M 2 )(1 - zf) (36) Gj+(*,t') - 2q 2 ^ T ±± / te a J 7 ^L-d( t _ 1 (z a )Pj.(4), 



2 7 s t'(t'-t) 

to identify the relevant kinematic dependencies. 1 

We start with the projection onto jtt —> jit partial f 2 2 J' + 1 f 1 — z s 7 

waves, which can be written as G ; j ,(s,t ) = 2q - J dz s ^ — ^i,-i( Zs ) 



-1 



00 

. 00 

/. 7 , + ( S ) = A+( S ) + -Jds'J2 Ktj-,(s,s')Imfj,, + (s') 

Ml J ' = 1 



x r*™ (38) 
1 z t 



00 Explicit expressions for J,J'<2 as well as the modifi- 

— / dt' G H — (s t')lmh t (t') cations for the subtracted case are given in App. [21 It 

7r J T ,„ JJ ' ' is important to note that while we have consistently sup- 

pressed isospin indices for the partial-wave amplitudes, all 
kernel functions are independent of isospin. 



J' even 



AMI 

00 



1 r 

/,,_(«) = Nj(a) + ~ J ds'J2 {Kj+(s,s')Imfj,, + (s') 

Ml J ' = 1 

+ Kjy,(s,s')Imfj,,.(s')) 

OO 

+ 1 j dt ' E (Gj+(s,t')lmhj,, + (t') 

AMI J ' CV0 " 

+ Gjy,(s,t')Imh JI ,. (*')), (37) 



3.4 £-channel projection 

Similarly, the projection onto 77 — > 7T7t amplitudes has 
the form 



1 /* 

J '° VCn = Nj(t) + ^jds'^2 (G++(t,s')Imfj,, + (s') 



Ml 



where N^(s) includes Born terms and — in case subtrac- 
tions were performed — pion polarizabilities. The kernel 1 



+ G+7,(M')Im f r ,_ ( s ')) 

OO 

I J dt ' E (Al+(i,t')Im^,+ (i') 



4 Note that and t' depend on the integration variable as 4M 2 J ' even 

well as on s and t. The subscript indicates evaluation at 

a = Ml and t = 0. ' + A+j" (t, t')Im hj',-(0 
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00 _^ 

hj,-(t) = NJ(t) + - I ds> Gj+(t,s')1mfj, i+ (s') 

Ml J ' = 1 

00 

+ ~ J dt ' Kjj;(t,^)Tmhj> t -(f), (39) 



J' even 



where, in the unsubtracted case, 

1 

g#'(M0 = J„, 2 2J o +1 / dz t (t'-t)Pj(z t ) 



8s'q' 2 2 



1 



1 



1 + z' s \ s' — s s' — u s' — a 

1 



2q' 2 2 

-1 



1 



1 



1 — zi I s' - s s' — u s' — a 
1 



+ o T' 1 1 

*#(*.*0 = 17777" / d 2t P^t)P/<(*0, 



- 1) 2 



JJ,[ ' j ~ t>(t>-4M%) 2 



dz t P/(z t ) 



1-zf 



Gj+(*y) = 



i(t - 4M 2 ) 2 J' + 1 
8s'q' 2 2 



d« t (l - zt )di (z t ) 



1 



KjJ,(t lt ') 



I + z' s \ s' — s s' — u s' — a J 

tft - 4M 2 ) 
t'(f'-4M"2)(f'-t) 

2^+1 K „ 2, ,, 



(40) 



d2 t (l - zt)di (z t ) 



i-z'r 



Explicit expressions for J, J' < 2 are provided in App. [B] 



3.5 Threshold and asymptotic behavior of the kernel 
functions 

In order to check our kernel functions and determine 
the convergence properties of the dispersive integrals, we 
study the behavior of the kernels at threshold and for 
s',t' ->■ 00. 

Based on (|2~T|) , one can show that for t —> 



H ++ {t,Ml) = h ,+ (t) - -h 2 , + (t) + (D(t 3 ), 
H + _(t,M*) = -f^ 2 h 2 At) + 0(t 2 ), (41) 



# subtractions 







1 




2 


K++(s,s'),K++(s,s') 


0(s'- 






-4\ 


C(s'- 5 ) 


K++{s,s'),K++(s,s') 


0(s'- 




0(s'- 


-4\ 


C(s'- 5 ) 


G+f(s,t') 


o(t'- 


- 3 ) 


o(t'- 


-4\ 


C(t'" 5 ) 


G&M 


o(f- 


- 4 ) 




-4\ 


0(t'~ 5 ) 


G 21 + (t, s'), G 22 + (t, s') 


0(s'- 


- 3 ) 


0(s'- 


-4\ 


0(s'- 5 ) 




o(t'- 


- 3 ) 


o(t'- 


-4\ 


0(t'~ 5 ) 



Table 1. Asymptotics of the kernel functions in the equations 
for /j,+ (s) and hj-(t). 



# subtractions 







1 




2 




K^+(s,s'),K^+(s,s') 


C(s'- 


- 2 ) 


C(s' 


- 3 ) 


G(s'- 


-4\ 


K^+(s,s'),K^ + (s,s') 


C(s'" 


- 4 ) 


0(s' 


- 3 ) 


0(s'- 


-4\ 


Krf(s,s'),K^"(s,s') 


C(s'- 


- 2 ) 


0(s' 


- 3 ) 


0(s'- 


-4\ 


K2f(s,s'),K^f(s,s') 


C(s'- 


- 3 ) 


0(s' 


- 3 ) 


0(s'- 


-4\ 


G w +(s,t'),G^ + (s,t') 


£>(*'" 


- 2 ) 


o(t' 


- 3 ) 


o(t'- 


-4\ 


G 2 o + (s 5 ^') 5 G 22 + (s, t') 


£>(*'" 


- 3 ) 


o(t' 


- 3 ) 


o(t'- 


-4\ 


G- 2 -(s,t') 


e>(*'- 


- 2 ) 


o(f 


- 3 ) 


o{f- 


-4\ 


G 2 2 ( s jO 






o(t' 


- 3 ) 


o{t'- 


-4\ 


Gqi (t, s'), Gq"^ (t, s') 


0(s'- 


- 2 ) 


0(s' 


- 3 ) 


0(s'- 


-4\ 


G21 {t> s ')i G%2 (t, s') 


0(s'- 


- 4 ) 


0(s' 


- 4 ) 


0{s'' 


-4\ 


A" + + (t,t'),^o + 2 + (*-*') 


o{f- 


- 2 ) 


o(f 


- 3 ) 


o{f- 


-4\ 


^ 2 + 2 + (*,f) 




- 4 ) 


o(t' 


- 4 ) 


o(t'- 


-4\ 





o(t'- 


- 2 ) 


o(t' 


- 3 ) 


0{t'~ 


-4\ 



Table 2. Asymptotics of the kernel functions in the equations 
for /j,-(s) and hj t +(t). 



where the hat indicates the subtraction of the Born terms. 
By comparing to ([21]) and (IB.ip . (jITj) implies that 



0{t 2 



K 



+±,(1) 
j j' 



5 G 
5 -K. 



-±,(2) 
2J' 



Go^ v ^(t, S ') 



T 2J' 
-^2J' 



(t,0 = O(i 2 ) 



2J' 

0(t 2 ) 
C(t 2 ^ 



+±, (2) 



(i,f) 



(t,0-O(t 2 ), 
0(t 3 ), 

G(t 3 ), 



2 ^ (2) (t,.s')=0(t 2 ), 
• 2 - J 7' (a) (t,0 = O(< a ), (42) 



where the superscript refers to the number of subtractions. 
We have checked that the explicit expressions in App. [B] 
fulfill these relations. 

Furthermore, is clear from (|37| and (l39|) that the de- 
pendence of the kernel functions on s and t must repro- 
duce the correct threshold behavior of the partial-wave 
amplitudes 

fj, + (s)=0(d 2J ), / J ,_(s)=0(q 2J ), 

hj,+(t) = 0(q 2 t (q tPt ) J ), hj^(t) = 0((q tPt ) J ). (43) 

The additional factor of q 2 in hj t +(t) is a manifestation 
of Low's theorem for low-energy QED |49j . which requires 
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Mi) 



Fig. 1. Box graphs constraining the boundaries of the double 
spectral functions. 



the full scattering amplitude to be equal to the Born terms 
at the threshold for Compton scattering. We have checked 
explicitly that the expressions provided in Apps. lAlandlBl 
indeed fulfill 

K++(s,8') = 0(q 2J ), G+J,( S ,t') = 0(q 2J ), 
KJ±( s , s ') = 0(cf J ), G,7±(M') = 0(q 2J ), (44) 
G+f,(t,s') = 0(qUqtPt) J ), K+f,{tJ) = 0(qUqtPt) J ), 
Gj+(t,s>) = 0{(q tPt ) J ), K-J,{t,t') = 0{(q tPt ) J ). 

Similarly, the asymptotic behavior of the kernel func- 
tions for s' — > 00 and t' — > 00, respectively, determines the 
convergence properties of the dispersion integrals. In par- 
ticular, one can directly read off which rate of convergence 
can be achieved when working with a certain number of 
subtractions. The corresponding behavior of the kernels 
for large values of the respective integration variable is 
summarized in Tables [T] and [21 Although in some cases 
the leading power vanishes, one can see that in general 
the kernels for /,/._(s) and hj, + {t) drop as the second, 
third, and fourth power in the integration variable, while 
the integrals related to /j, + (s) and hj,^.(t) will converge 
one order faster. 



4 Domain of validity 

The original Roy equations for tttt scattering have been 
shown to be rigorously valid up to s max = 60M^ based 
on axiomatic field theory [2J. This range can be extended 
further to s max = 68M^ if the nn amplitude T fulfills 
Mandelstam analyticity |50j . i.e. it can be written in terms 
of double spectral functions p 8U , pt u , and p s t as 

i-r-f ,\ 1 f f A 1 A 1 Psu{s',u') 

r (M) = ZS // dsdu 7^" 



1 



(V — s)(u' — u) 
dt'du' ^''^ 



(t'-t)(u'-u) 



(s'-s)(t'-t) 



(45) 



In this section, we derive the corresponding limits on our 
system of jtt — > 77r Roy-Steiner equations. For more de- 
tails of the derivation we refer to similar work on ttK and 
nN scattering [46ll48l[5T] . 

The central objects of the discussion will be the bound- 
aries of the support of the double spectral functions, which 



determine the integration range in (j4"5"]) . These boundaries 
can be inferred from the box diagrams depicted in Fig. [T] 
These diagrams are to be understood as generalizations 
of four-propagator box diagrams |52j . with one or more 
lines replaced by a particle whose mass is equal to the in- 
put mass of the lowest-lying intermediate state accessible 
to the interacting particles. If we neglect the possibility 
of photons in intermediate states, the pertinent states in- 
volve one and two pions. The diagram (I) represents the 
mechanism that produces the boundary of the support 
of p s t, while (II) and (III) are relevant for p su . "Double- 
spectral regions" are then defined as the portions of the 
Mandelstam plane that obey s + t + u = 2M% and where 
any of the functions p s t, p su , Ptu has support. 
From (I) we find that 



h{s,t) - t(s - mil) - 4M%(s - Ml) = 



(46) 



defines the st boundary of support, and the same with 
shii for the tu boundary, while 

MM) = &m(M) = su + Ml(9t - Ml) = (47) 



is the boundary of p s ,, 
may be described by 



In particular, the boundary of p s t 



t = T st (s) = T I {s), 



(48) 



where t = T\(s) follows from solving b\ (s, t) = for t, and 
similarly for p tu and p su . 

These domains of support restrict the range of validity 
of the Roy-Steiner equations in two ways: 

1. The partial- wave expansion of the imaginary parts in 
the dispersion integrals converges only for scattering 
angles z that lie within the large Lehmann ellipse [5.3, ■ 
This ellipse can be constructed as the largest ellipse 
in the complex z-plane which does not reach into the 
double spectral regions. Given a value of a, this con- 
straint can be translated into an allowed range for the 
parameter b, since ([55)1 relates b to the angle z. 

2. A specific value of b is only allowed if the hyperbola 
(s — a)(u — a) = b does not enter the double spectral 
regions. 

The allowed values of b (for a given a) must respect both 
constraints and must do so in the integrals which occur in 
expressions for both the 771 — ¥ "/it and 77 — ¥ tttt ampli- 
tudes. Once we have a constraint on the allowed values of 
b we can translate it into limits of the domain of validity 
of the full Roy-Steiner system because of the partial- wave 
projection of the dispersion relations. For example, in the 
s-channel, we need — 1 < z s < 1, such that the maximally 
allowed s max is the largest value of s for which both 



and 



(s-a)(2Ml 



(s - a) 2Ml -s-a + 



a) 



{s-Mlf 



(49) 



(50) 



lie within the allowed range for b. In a similar fashion 
the requirement that < z 2 < 1 determines i max in the 
i-channel projection. 
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We begin with the Lehmann-ellipse constraint in the 
s-channel. The partial-wave expansion converges in an el- 
lipse with foci at z' s = ±1 



(Re^) 2 , (Im^) 2 



-4? 



= 1, 



(51) 



so that the semimajor and semiminor axis A s and B s are 
related by 

Al-B 2 s =l. (52) 

The maximal value of z' s that does not enter the region 
where p st ^ is given by 



1 



2s'T st (s') 
(s' - M 2) 2 ' 



(53) 



such that the st boundary provides the Lehmann-ellipse 
constraint 

- zT x <z's< zT x - (54) 
Translating this into a restriction on t' we find 



T' at {s')<t'<T st {s'), 

T'As') = - (s/ ~f' )2 -T st (s'), Vs' G [A/ 2 , TO ), (55) 



in the s-channel integral of the dispersion relations. As 
(j33|) defines a linear relation between b and t' , (l55l) trans- 
lates into a condition on b 

b:(s',a)<b<b+(s',a), (56) 

where 

6 s +(s', a) = (*' - a)(2Af 2 - s' T' st {s') a), 
b-{s',a) = (s' - a) (2 A/ 2 - a' - T st (s') - a). (57) 

We may then define 

(58) 



b s (a) = max b s (s' , a) 

s'£[MJ,oo) 



to be the minimum/maximum of bf(s',a) within the do- 
main of integration. Similar restrictions are provided by 
ptu and Psu, and the intersection of the resulting con- 
straints on b defines the limitations imposed by condition 1 
due to the s-channel part of the dispersion relation. 

The same condition in the case of the ^-channel reac- 
tion involves a slightly more complicated argument, be- 
cause in this case the relation between the CMS angle z' t 
and b is not linear 



, 2 _ (f - 2M 2 + 2a) 2 



46 



(59) 



Consequently, we need to consider the Lehmann ellipse for 
(Rezf-I) 2 ; (Imzf " 



-/ 2 



Af 



+ 



Bf 



= 1, 



(60) 



where the parameters are related to those of the ellipse 
for z' t by 

A t = A\- X -, B t = A t B t . (61) 

Rewriting (|46[) in terms of v and t and inserting the result 
in (IT5j) , we obtain the boundary of the double-spectral 
region in terms of z' t 



N(t>) 



Hv't N{t ' ] t'-AM* ' 
such that p s t imposes the restriction 

bt(t',a) <b<b+(t',a), Vt' € [4M 2 ,oo), (63) 

with 

bt{t',a) = \{t' -2Ml + 2af-\N{t')\ 

b+(t>, a) = - 2Af 2 + 2a) 2 + ±N(t' f - 4q^ 2 , (64) 

and similarly for p tu and p su . Defining the variables b^(a) 



t'(t' + 12Af 2 ) 



(62) 



bt(a)= min bt(t',a), 

t'G[4M2,oo) 

K(a) = , max bt(t',a), 



(65) 



the constraints 

b-{a)<b<bt{a), bt(a)<b<b+(a) (66) 

then provide, for a given a, the range of allowed values for 
b that satisfies the Lehmann-ellipse constraint for both the 
s-channel and t-channel integrals. 

If the associated hyperbolae in the external variables s 
and it, (s — a)(u — a) = b, do not cross the double spectral 
regions either, then this determines the kinematic regime 
in which the partial-wave projection is valid. 



4.1 s-channel projection 

The freedom in the choice of the parameter a in the con- 
struction of the hyperbolic dispersion relations can be used 
in order to maximize the region in s within which these 
dispersion relations are valid. 

As a first step, we consider the restrictions due to p st in 
the s-channel part only. The strategy to find the optimal 
value of a goes as follows: the Lehmann-ellipse constraint 
imposes the condition that all allowed values of b must 
fulfill b~(a) < b < bf(a). As t varies only within 



( s ~ Af 2 ) 2 



< t < 0, 



(67) 



this limits the range of values of b that are needed for a 
given s. The value of s which is the maximum one possible, 
Smax, will be that which yields a result where the smallest 
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b coincides with b s (a), and the largest b coincides with 
bf(a), i.e. the value of s such that the solutions of 

(s-o)(2A^-»-o) = 67(o), 

<»- M *> 2 W«0 (68) 



(s - a) 2M; - s - a 



coincide. This procedure results in 



a = -41.3 M; 
6+ (a) = 2852 M*, 



s majc = 27.8M^, 
6-(a) = 1071M* 



(69) 



/9 SU and pt u do not yield additional constraints. 

The investigation of the i-channel contributions pro- 
ceeds along the same lines: demanding that the solutions 
of the two equations 



{s-a)(2Ml 

■2 .. , ( s " 



(s - a) \2M* - s - a 
coincide, we find 



- s — a) 
M 2 ) 2 



= b t ( a ), 

= b+(a) (70) 



I 2 

1 TV > 



bt(a) = 308.4 Af 4, 6" (a) = —298.1 M*. 



(71) 



Thus, the s-channel constraint (j6"9")l is weaker than the 
i-channel restriction ([7T|) . as can be deduced from Fig. [21 
where the situation for a = —9.8 M 2 is displayed: both 6/" 
and 6^~ also lie within the allowed region for the s-channel. 
Since the hyperbolae resulting from the choice (fTTj) do not 
enter the double spectral regions either (see Fig. [3]), ([7T|) 
constitutes the final answer for the range of validity of the 
s-channel projection. 



1500 



1000 - 



500 



-500 



-1000 - 



-1500 L 





V 6+(s', a) ^— - 











1500r 



1000 - 



500 - 



-500 



-1000 - 



-1500 



10 



15 



20 



25 



30 



35 



40 



Wo) 




20 25 

f I M% 

Fig. 2. Allowed range of b for a — —9.8 M% in the s-channel 
projection. The red lines correspond to bj~ (a) = 308.4 and 
67(a) = -298.1 Mi, respectively. 



4.2 i-channel projection 

As we may start from a different set of hyperbolic dis- 
persion relations to perform the t-channel projection (|3"9")l . 
there is no need to use the same value of o as in the s- 
channel. To perform the t-channel projection, we need to 
consider 



< zt = 



(t - 2M 2 + 2a) 2 - 46 



16q 2 p 2 



< 1. 



(72) 



Similarly to the s-channel case, we determine the optimal 
choice of a such that the solutions for t with z? = and 



1 are the same. This yields 



a = -7.5 Ml, 
b+(a) = 298.4 M* 



-17.4 M 2 < t < 51.6 M 2 
67(a) = -316.8 M*. 



lGeV 2 , 
(73) 



Again, this region is driven by the constraint on the t- 
channcl Lehmann ellipse, which provides a stronger con- 
straint than that involving b^(a) and bj(a). We will use 
the value a = —7.5M 2 in the rest of the paper. 

As we eventually aim to investigate the properties of 
the a, we also need to consider the domain of convergence 



40 






20 


Ptu \J \ 


\ Psu 







3 






-20 




Pst 


-40 








-40 -20 


20 40 




s 1 Ml 




Fig. 


3. Double spectral regions and 


hyperbolae for a = 


-9.8 Ml, b+{a) = 308.4 and 6 t _ (o) = 


= -298.1 M*. 
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■20 



-10 



10 20 

Ret / A/, 2 



30 



40 



50 



Fig. 4. Domain of validity in the complex t-plane for a — 
— 7.5 M%. The blue dots refer to the position of the cr-pole 
at t = (6.2 ± il2.3) Ml and the red dots to the f at t = 
(51.4±il.4)M^ p. 



in the complex plane. We will restrict our analysis to the 
value of the hyperbola parameter a = — 7.5M 2 . Then, the 
constraints from both s-channel and t-channel Lehmann 
ellipses yield ellipses of allowed values for b in the com- 
plex 6-plane. As the ellipse for the s-channel contains the 
ellipse for the i-channel, it suffices to consider the latter. 
The resultant region of the complex t-plane is depicted in 
Fig. |U It safely encompasses the position of the a pole, 
and marginally that of the /o(980). 

The result that the equations for the t-channel are rig- 
orously valid up to i max = 1 GeV 2 seems to shed doubt on 
recent dispersive fits [30) to high-statistics 77 —¥ irir data 
[8 , 9J in the energy region up to 1.28 GeV, where a lot of 
effort was put into building a good description of K K dy- 
namics above 1 GeV. One possible way to obtain a higher 
upper limit on the range of validity of the 77r Roy equa- 
tions would be to relax the assumptions about the bound- 
aries of the double spectral functions: if one assumed that 
the spectral strength of the 2-pion intermediate states de- 
picted in Fig. [T] only set in at an effective mass larger 
than 2Mn, the domain of validity would be extended ac- 
cordingly. For instance, a threshold of m c s = SM^ would 
produce t max = 82.8 M 2 = (1.27 GeV) 2 at a = —12.0 M 2 . 



5 Muskhelishvili-Omnes problem for 

77 — > 7T7T 

We now turn to the resolution of the equations for 77 — > 
7T7T. We truncate the system at J — 2 both for s- and 
i-channel contributions. The generalization to higher par- 
tial waves is straightforward. In this approximation, the 
equations can be brought into the form 



i(0 



. , t 1+i [ 

(t) + —J 



dt 



4A/? 



, Imfo 0;+ ft') 



h 2 .+{t) = A 2 .+{t) 

t 2 (t-4Ml) 



dt'- 



Imft 2 ,+ (*') 



,7 j t' 2 (t'-4M2)(t'-t)' 

AMI 



& a ,_(t) = 4;>_(t) 

t 1+l (t - 4M 2 ) 



(74) 



dt'- 



Im/i 2 ,_(i') 



j t n + l (t> -4M2)(# -t)' 

4M2 



where i £ {0,1,2} indicates the number of subtractions, 
and Aj t ±(t) includes the Born terms, subtraction con- 
stants, and integrals involving the imaginary part of other 
partial waves 



4: + W^o +W (i) + ^ j dt^K^^>(tX)lmh 2 Mt') 

OO 

Ml J'=1.2 

+ G+7 ( %/)Im/ J ,,_( S ')), 

OO 

A 2 , + (t) = N+(t) + - f ds' ]T (Gi+{t,s')lmfj, t+ {s') 

+ G+7(t, S ')Im/j',-(s')), 
4?-(*) = # a ~ (<) (*) 

OO 

+ lJ ds 'H G- J , W (t, S ')In 1 /.7', + ( S ')- (75) 

Several comments are in order: first, the equation for 
h 2i+ (t) is not affected by subtractions. Second, the equa- 
tions for the D- waves decouple, as the corresponding ker- 
nels relating them to the S-wave or to each other van- 
ish. Conversely, the S'-wave does not decouple from the 
D-waves: both D-waves are needed as input to solve for 
/io,+(i). The consequence of these observations for the nu- 
merical evaluation is that we first solve the equations for 
the -D-waves separately, and then we use these solutions 
to compute as input for the S'-wave. 

Assuming elastic unitarity 

Jmh Jt ±(t) = <r{t)h l± (t)tj{t)* (76) 
with 7T7T partial waves tj(t) in the normalization 

e 2i5.,(t) _ 1 



-(i), 



1 



tj(t) 



2ia(t) ' 



(77) 



the phase of hj t ±(t) coincides with the 7T7t phase Sj(t) 
("Watson's theorem" [54]), and ([74]) takes the form of a 
Muskhclishvili-Omnes problem [T8"UT5] . whose resolution 
we will review in the following section. 
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5.1 Muskhelishvili-Omnes problem with finite 
matching point 

As we have argued in Sect. HJ the domain of validity of 
(j74|) is restricted to the energy range below 1 GeV. We will 
follow here the strategy of , namely to assume that 

the imaginary parts of the amplitudes are known above a 
matching point i m , and solve the equations in the regime 
between threshold and t m . The input that we will use both 
for the s-channel contributions and the high-energy regime 
above t m will be discussed in detail in Sect. [7] 

In this finite-matching-point setup, the solution in 
terms of Omnes functions reads 

s OO 

ho,+{t) 4>,+(*) + — — J & t ,i + i {t >-t)\n (t>)\ 



, smS (t')A^ + (t') 



dt 



iMz 



t'i+*(t'-t)\n (t')\ ' 



ha,+(t) = A 2 ,+(t) + 



t 2 (t-4M*)n 2 (t) 



x( dt' 



4MB 



sm5 2 (t')A 2 , + (t') 

t> 2 (t>-4M%)(t>-t)\n 2 (t>)\ 



dt' 



t>2(t>-4MZ)(t'-t)\f2 2 (t>)\ [' 



h 2A t) = A 2 Ut) + tl+t{t - 4M2Mit) 



xl dt' 



sin S 2 (t')A^_(t f ) 



v LU: 



dt' 



where 



t>i+i(t> -4M%)(t> -t)\n 2 (t>)\ 



Imh 2 ,-(t') 1 

t>i+i(t>-4M%)(t>-t)\n 2 (t>)\ [' 



^ 4MB 



t) 



\ko,+ [t)\ = cos<5 (t) + 



i 1+ i^oWI 



x <^ + dt 



t'i+*(t'-t)|0,(t')l 



'4MB 
00 



/■ Imfo , + (t') I 

J t'^(t'-t)\Q Q (t')\V 



(78) 



(79) 



Defining the physical amplitude by the limit t — > t + ie, 
(|78|) can be rewritten at@ 



(80) 



The symbol f indicates the principal value of the integral. 



and similarly for the D- waves. There are several subtleties 
related to the behavior of |J?j(t)| for t — > t m 48 



I'M*) I 



t„ 



t 



, ij('m) 



(81) 



in particular the integrals in ([50)1 diverge for Sj(t m ) > ir. 
In such a case, there are non-trivial homogeneous solu- 
tions whose coefficients can be used to absorb these diver- 
gences, but, of course, the presence of such solutions in- 
troduces undetermined constants in our result for f2j(t). 
For instance, the solution for h 0:+ (t) for n < S (t m ) < 2tt 
involves one free parameter a, 



iv+wi = 4?+ (*)«*$>(*) + 



t i+i \n (t)\ 



a + t j- dl 

4A/B 



+ t dt 



(t m - t)lt 

,(t m -1/)sm6 Q (tf)A^ + (f) 
t>*+i(t>-t)\Q (t>)\ 

t')lmh ,+(t')} 



t»+i(t>-t)\n (t>)\ 



(82) 



a can be fixed if one assumes, in addition, knowledge of 
the derivative of ho t +(t) at t m . Conversely, the fact that 
the phases of the 1 = 2 partial waves are negative induces 
different complications, which we will address in Sect. 15.21 
The part of the integrals in ((78)) involving the pion 
polarizabilities can be explicitly performed based on the 
spectral representation of the Omnes functions 



nj l (t) 

Q-/(t) 



dt' 



IMz 



1 - - 



{t'-t)\Qj(t')V 
sm5j(t') 



dt' 



iM'i 



Qj 1 ^) = 1 -tf2j(0) 



t'(t'-t)\Qj(t')Y 

sin 5 j ft') 



(83) 



t l 

7T 



dt' 



4MB 



t^(t'-t)\Qj(t')[ 



where the dot denotes the derivative with respect to t. 
The results of this modification, which, in particular, only 
involve integrals over 



Aj !± (t) = A J<± (t)-ANf(t), 



(84) 



are summarized in App. IC.ll (with ANf(t) defined in 
App.l]). 



5.2 Sum rules for I — 2 

If 5j(t m ) < 0, f?o(t) diverges at t m , such that (|78l) breaks 
down. However, we may rewrite the solution as 

v +(t ) = 4",w + tW "° ( ' )(t -- t) 
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^(tm -m> -t)\n (t')\ 

lmh Q ,+ (t') \ 

t> l +*(t m ~t>)(t< -t)\n (t')\j 




(85) 



the uniqueness properties of the solutions of tttt Roy equa- 
tions presented in 55,56,57. There, an additional free pa- 
rameter occurs each time the phase at the matching point 
crosses an integer multiple of tt/2. Indeed, to derive this 
result the Roy equations are linearized, and in the one- 
channel approximation one finds that the difference 



x _ S'js) - 6(3) 
a(s) 



(89) 



where 



between two solutions S(s) and S'(s) for the tttt phase shift 
must fulfill [56] 



t i t*-"p-t m )\n (t>)\ 

AMI 
00 

tl +i f dt , Imfto,+(f) 



TT J t'i+*{t> -t m )\f2 (t>) 



(86) 



Demanding continuity at the matching point implies that 
— 0, since otherwise ho,+(t) would diverge at t m . In- 
deed, multiplying the first equation of ([75)1 with {Iq 1 ^) 
and evaluating it at t = t m explicitly proves that cS 1 ^ = 
must hold in order for /io,+(i m ) to be finite. The final re- 
sults for all amplitudes as well as the explicit form of all of 
the possible sum rules which could be obtained are given 
in App. El 

One case that is of particular interest is the 1 = 2 
S-waves, where the tttt phase shift is substantial, and neg- 
ative, at the matching point. Thus, we obtain from the 
once- and twice-subtracted versions of the Muskhelishvili- 
Omnes representation 

0=^< m (a 1 -/3 1 ) / = 2 + /( 1 ), 
la 

lot 

+ 2^ (a2 ^ 2)/=2 + /(2) ' 



(87) 



where 




dt' 



lmh 0>+ (t') 



t ,1+i (t'-t m )\n (t')\ 



ie {1,2}. 



As we shall see below, the second sum rule, in particular, 
provides a useful and novel constraint on (0*2 — /?2) /=2 - 



5.3 Uniqueness and comparison to rnv Roy equations 

The pattern in which free parameters emerge in the 
Omnes solutions is reminiscent of the results concerning 



cos25(s)(f>(s) = 



4Mi 



ds' 



4 A/2 



sin2<5(s')<M s ') 
(s' - 4M%)(s' - s)' 



(90) 

where s m is the matching point. In the tttt case the phase 
shift is the quantity that is to be determined, while in 
77 — > tttt that phase is input and the modulus of the 
amplitude \h(t)\ unknown, but, as we shall now see, the 
mathematical structure is the same up to a factor of 2. 

As an example we consider the unsubtracted equation 
for the 5-wave. In 77 — ¥ tttt the difference 



between two solutions of ([71)1 obeys 



cos<?o(t)V>o,+ (*) = 



dt 



t'(t' - t) 



(91) 



(92) 



The presence of sin5o(t')) rather than the sin2<5(s') of the 
tttt case, explains why the multiplicity of solutions in tttt 
scattering is |_2J(s m ) /7rJ instead of |_<5o(*m)/7rJ in 77 — > 
tttt, where \x\ denotes the largest integer < x. 

In addition, there is also an analog of the sum rules 
discussed in Sect. 15.21 in tttt Roy equations. In the single- 
channel approximation <fi(s) vanishes if S(s m ) < tt/2. If 
the effects of the coupling to other partial waves are taken 
into account, the phase-shift difference for a channel i with 
Si(s m ) < tt/2 may be written as |57j 



<j>i(s) = ( S -4M 2 )G I (s) J E I (s), 



Gi{s) = exp < ^ J ds'^ 

^ AMI 



(93) 



where Hi(s) includes information on other partial waves. 
Once again, Gi(s m ) diverges if Si(s m ) < 0, such that 



Hi(s m ) = 



(94) 



in this case. The difference as compared to 77 — > tttt is that 
in tttt the Omnes representation is only available for phase- 
shift differences and not for the phase shifts themselves. 
Thus, the constraint manifests itself rather subtly by re- 
ducing the number of free parameters in the manifold of 
solutions. However, the mathematical input that leads to 
this constraint, namely continuity of the Muskhelishvili- 
Omnes representation at the matching point, is the same. 
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5.4 Comparison to previous work 

The key result of our derivation of dispersion relations for 
77r scattering and 77 — > irir based on Roy-Steiner equa- 
tions is that there is a term that couples S- and Z?-waves, 
cf. (175[) . In all previous dispersive treatments of 77 — > irir 
each partial wave was considered separately. Moreover, 
to the best of our knowledge, this is the first time that 
a finite-matching-point representation has been employed 
for 77 — > 7T7T. The practical consequences of both these de- 
velopments will be discussed in Sects. [7] and \E\ However, 
we can consider the limit t m — > 00 in our equations in 
order to delineate the differences between our formalism 
and the recent works [TT)ll2"Tll3Tj] . which also used disper- 
sion relations to attack the problem 77 — > irir. 

In [TS] , a once-subtractec0 dispersion relation for the 
S'-wave is considered. The subtraction constant is fixed 
by assuming ho,+ (t) oc to(t) and using ChPT information 
on the Adler zero of the 7T7t amplitude. This representa- 
tion depends quite strongly on the details of the 7T7t phase 
above the KK threshold already at energies > 0.5 GcV. 
For this reason, another subtraction was performed in [27] 
at the energy t\ where the irir phase crosses ir, the sub- 
traction constant being fixed by the requirement that the 
cross section at t\ does not become outrageously large. 
The Born terms as well as vector and axial-vector reso- 
nances were used to approximate the left-hand cut. 

In [3D], an Omnes representation for S- and D-waves 
is constructed that explicitly takes into account the KK 
channel, and in addition includes tensor resonances in the 
description of the left-hand cut. For the S- waves two sub- 
tractions are performed, while the D-waves are treated 
differently for the two isospin channels: for I = 2, no 
subtraction constants for /i2,+ and /12,- are provided (in 
our conventions this corresponds to the unsubtracted case 
for /12,-), while for I = an additional subtraction in 
both partial waves is performed. Thus the treatment of 
1 = corresponds to our once-subtracted case for /i2,-. 
But the Roy-Steiner analysis shows that the subtraction 
constant for li2,+ cannot be related to dipole and quadru- 
pole polarizabilities, as the equations for this partial wave 
are not affected by the corresponding subtractions. For 
this reason, the subtraction constants in [30] determined 
by fits to data (together with several chiral constraints) 
can be translated into pion polarizabilities, but in general 
the converse is not possible. We have checked explicitly 
that our results agree with |30j in the limit i m — > 00, once 
the KK channel is switched off and the additional sub- 
traction in /i2,+ dropped. 



6 Photon coupling of the a resonance 

We define the airir coupling constant such that the 
full isospin 1 = irir scattering amplitude on the second 

6 Note that in this work we do not count powers of t' or 
t'— 4M% that are present for kinematical reasons alone and thus 
do not require any subtraction constants. This is not always 
the convention employed in the literature. 



Riemann sheet Tj°j near the position of the a pole 

t a = (u a - i^j (95) 

can be written as 

00 n 2 
T° = 327r$>J+l)t° II (i)P / ( Zt ) = fzz- (96) 

As the cr-pole occurs in the S-wave, all other partial waves 
only contribute to the background. Similarly, in 77 — > irir 
we take 

2 

e H ++,u = — t _ t — ■ (97) 

In these conventions, the widths in the narrow-width ap- 
proximation are 

I 2 /, 4M| _ 7ra 2 |g J77 | 2 

CT ™ ~ 32irM a V Ml ' CT77 ~ M a ' 

Note that since the large strong width of the a renders 
the applicability of these relations questionable, we use 
<7 CT7r7 r as determined from the residue of the pole. Although 
the formula for -T CT77 and <7 CT77 suffers from the same de- 
ficiency, it is conventionally employed in the literature to 
illustrate the relation between the two quantities. As the 
direct determination of T ail from the position of the pole 
is not possible in view of the large strong width, we follow 
this convention. The coupling constant itself can always 
be recovered by means of 

The analytic continuation of h® , (t) into the complex 
plane is given by [27J3 

h°o, +t iM = (1 - 2^(^,„(*))< +iI (i), (99) 

such that 

k---{ 2 M) 2{h ^ < ioo) 

where h® , is evaluated on the first Riemann sheet. As- 
suming the position of the a pole and its coupling constant 
to two pions is known, we can thus infer 3o- 77 (and hence 
_r, T77 ) from the value of the 1 = S'-wave of 77 —¥ irir 
evaluated at t a . 

7 Input 

To solve the Roy-Steiner equations for the 77 — > irir par- 
tial waves we must specify the input for Im/j i= |-(s) in the 

7 We neglect 77 intermediate states in the unitarity rela- 
tion. This is the same approximation as used in |T], where the 
a pole is deduced from irir scattering with electromagnetic in- 
teractions switched off. Since the width corresponding to the 
77 channel amounts only to a few keV, the ensuing shift of 
the pole is much smaller than the uncertainty of its position 
as quoted in [T|. 
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whole energy range and for Im hj t ±(t) above the matching 
point. One could, in the spirit of the Roy-Steiner analysis 
for ttK scattering [38], consider the equations for s- and 
t-channel partial waves simultaneously, and determine a 
solution of the whole system by an iterative procedure. 
However, the Roy-Steiner equations for 77r — > 777 are less 
powerful than those for ttK — > ttK, since the only con- 
tributions that can be obtained in the s-channel without 
additional input are determined by elastic unitarity. In the 
77r case, 77r intermediate states are suppressed by e 2 and 
thus expected to be numerically negligible. By compari- 
son, in the ttK case, ttK intermediate states dominate the 
unitarity relation at least up to 1 GeV [48] . 

For this reason, we will drop jtt intermediate states 
altogether and content ourselves with the resonance de- 
scription of the imaginary parts of the Compton-scattering 
amplitudes constructed in [3D], where the resonance con- 
tributions are eventually included in a spectral representa- 
tion with an integration cutoff of —5 GeV 2 . In our frame- 
work the effect of the resonance description of Imfj t ±(s) 
on hj t ±(t) can be directly read off from (l39l) or ([75]) . 
cf. Sect. 17.11 Physically, one can understand this summa- 
tion of resonances as an effective description of multi-pion 
states in the s-channel for 777 — > 771", or, equivalently, in 
the i-channel for 77 — > tttt, amounting to approximating 
the multi-pion cuts by a sum of poles. This model could be 
improved upon at low energies by explicitly incorporating 
the 2- and 3-pion intermediate states, and using ChPT to 
constrain their interactions. However, these intermediate 
states enter the 777 amplitude at 0(e 2 p 6 ) and 0(e 2 p 4 ) in 
the chiral counting, respectively, such that we will leave 
their incorporation to future work. 

We now turn to the input for lmhj,±(t). We choose 
the matching point as 

\ft^ = 0.98 GeV, (101) 

which, on the one hand, ensures that Sa(t m ) < tt, avoiding 
a free parameter in the Omnes solution, and, on the other 
hand, extends the energy range as far as possible. As the 
cross section above 1 GeV is dominated by the /2(1270) 
resonance, we put hj^±(t) = above t m for all partial 
waves except for h^^it), which we match to a Breit- 
Wigner description of the /2(1270), cf. Sect. 17.21 As we 
will show in Sect. [8J this approximation allows for a rea- 
sonable description of the cross section. A more detailed 
description of this region is not necessary, because our 
subtracted dispersion relations emphasize the low-energy 
region which is the focus of this work. 

As input for the tttt phases we use the results of an ex- 
tended Roy-equation analysis of ttit scattering [58], which 
in particular ensures that the phases and the pole posi- 
tion of the a are consistent, since [S8J coincides perfectly 
with the older analysis [5DJ at low energies. (The impact 
of the high-energy region in tttt scattering on the cr-pole 
position was shown to be negligible in [JJ.) To estimate 
the uncertainties due to the tttt input, we also consider 
the tttt phases determined in a recent study of Roy-like 
equations [BDJ. The parameters of the a resonances cor- 
responding to both approaches, which we will refer to as 







r CT 


Scttttt/ V2 


CCL [T][6"T] 


441±* 6 MeV 


544± 4 ® MeV 


3.3 GeV 


GKPRY [62] 


4571^ MeV 


558±^ MeV 


3.59±g;^ GeV 



Table 3. Mass, width, and tttt coupling constant of the a. 





ChPT 38,39 


GMM [30] 




-1.9 ±0.2 


-1.25 ±0.17 


{ax 


1.1 ±0.3 


1.22 ±0.12 


(a a -I3 2 r° 


37.6 ± 3.3 


32.1 ± 2.1 


(a 2 + P2Y 


0.037 ± 0.003 


-0.19 ±0.02 


{ax-PxT* 


5.7 ±1.0 


4.7 


{ax+pxY* 


0.16 [0.16] 


0.19 ±0.09 


(a 2 - P2Y ± 


16.2 [21.6] 


14.7 ± 2.1 


(a 2 + fcY* 


-0.001 [-0.001] 


0.11 ±0.03 



Table 4. Dipole and quadrupole pion polarizabilities in units 
of 10 _4 fm 3 and 10 _4 fm 5 , respectively. The numbers in brackets 
refer to the LECs from [63] . 



CCL and GKPRY, respectively, are given in Table [3J and 
are consistent within error bars. 

Finally, our results depend on the input chosen for the 
pion polarizabilities. For definiteness, we will consider the 
two sets of parameters compiled in Table [3J which we refer 
to as ChPT and GMM, respectively. The polarizabilities 
in the isospin basis follow from ([2"D]) . Note that in [3DJ the 
dipole polarizability of the charged pion was only allowed 
to vary within the range of the ChPT prediction. Unfortu- 
nately, the charged-pion quadrupole polarizability c*2 — /?2 
is rather sensitive to low-energy constants (LECs): the 
first number, 16.2- 10 -4 fin , corresponds to the resonance- 
saturation model of [3DJ , while taking the LECs from [63 
yields 21.6 • 10~ 4 fm 5 . 



7.1 Resonances in 77r — >■ jtv 

We use the resonance model constructed in [3D], with 
the contribution of vector (V), axial- vector (A), tensor 
(T), and axial-tensor (Ta) resonances to the Compton- 
scattering partial waves in the narrow-width approxima- 
tion, to define the imaginary part of the 777 — ¥ "/tt ampli- 
tudes 

Im/X±(s) = ±^TrCy(m 2 v - M 2 ) 2 S(s - m 2 v )S JU 

Im/£ ± (s) = ^TrC A (m 2 A - M 2 ) 2 S(s - m 2 A )S, ni 

Im/J ± (,s) = ±Kc T {m ^~ 2 M ^ 6(s - m 2 T )5 J2 , 

t f T A / \ 2 (m T — M 2 ) 4 2 
lmfj±(s) = -ttC Ta — 2 5{s - to Ta )<5j 2 , 

(102) 
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where nij, i € {V, A, T, Ta}, denotes the mass of the res- 
onance, and the coupling constants Cj are related to the 
widths by 



3m v 



(m 



1- M l 
3m A 



i X — Ct(->T ; — E , Jt a =<»OTa 

5m T 



K A -^) 5 



5m T 



(103) 



Inserting (fT02|) into d39j yields 
2 



= 3^v(™ v - M*) 2 (£++(*, m v ) - G+r (t, m v )) , 



3 



#,+(*) = jjG A (m| - M 2 ) 2 (G++(t, ml) + G+ n (t, m|)) , 



G A (m A 



MiYG-^(t,mil 



(m T - M 2 ) 4 



hl_(t) = -Or 



(m T — -M 2 ) 4 ~ 



G++(t,m 2 T )-G+-(t,m 2 T )), 



t , s 2 (m Ti - M 2 ) 



2 Gj-ttt' m T); 



2^4 



Ta 



x (G+ 2 + (t,m| A ) + G+-(i,m| A ) 



5 C t a ^5 G J2 +(<,m TA ). 



(104) 



We include all resonances listed in [30]. We have checked 
that (|104p agrees with the results quoted in [3U] : once the 
ambiguous term linear in t in ho^(t) in j30] is removed, 
we recover that result by taking the limit a — > 00 of our 
unsubtracted kernel functions. 



7.2 Including the / 2 (1270) 

To incorporate the D-wave resonance /2(1270) we follow 
[24] . Starting from 

Ctpp = C^T^d^PdyP, £x 77 = e 2 C^ v F m F v a 

(105) 

to describe the coupling of a tensor resonance to pseu- 
doscalars and photons, respectively, we find 



fin r~f\ 



D = 



6(t — m^) 

K_y r-p V_ ' rp 

4(t-m|) ! 



A = -^r 4M? 4-4r -t S 



2^ 

4 



(106) 



such that 

#4 



= c^tW (107) 

4 t — TTIt 





full 


a — > 00 


no resonances 


CCL 


-0.62 


-1.15 


0.61 


J (1) , GKPRY 


-0.63 


-1.17 


0.60 


7< 2 >, CCL 


3.45 


3.58 


2.08 


7 (2) , GKPRY 


3.40 


3.53 


2.03 



Table 5. Integrals in the 1 = 2 sum rule for the full left-hand 
cut, in the limit a — > 00, and with resonances switched off. 



ai - Pi 



a.2 — P2 



total 



ChPT 1.03 ±0.14 
GMM 0.80 ±0.14 



-4.29 ±0.78 0.18 ±0.85 
-3.49 ± 0.60 0.76 ± 0.68 



Table 6. Individual contribution to (|87[) from the polari- 
zabilities (first two columns) and total value of the right-hand 
side of the sum rule (third column) . 



We see that a resonant contribution only occurs in h,2,- (t), 
while the non- resonant background in ho, + (t) can be dis- 
carded. Taking the full width of the /2(1270) into account 
and dropping the non-resonant background, we obtain 



4/2 

mJ 2 <J(mj 



5y/6 t - mj 2 + imf 2 r f2 



(108) 



In fact, in Sect. 18.21 we will restore the background in or- 
der to describe the cross section for 77 — > tt + tt~ above 
the matching point. The coupling constants can be deter- 
mined from the partial widths 



r 



(q- ) (m 



4Mi 



960tt 



7-1 2 / s~i*i \ ^ 3 



(109) 



For the f% parameters we use as input [64] 

m h = 1275.1 MeV, r h = 185.1 MeV, 

= 156.9 MeV, i/ 2 _>. 77 = 3.03 keV, (110) 

such that 

|G£| = 16.06 GeV~\ \C} 2 1 = 0.21 GcV" 1 . (Ill) 

However, the relative sign of the couplings cannot be in- 
ferred and must be fitted to experiment. 

8 Numerical results 
8.1 Sum rules 

Before examining the reactions of interest, we turn to the 
numerical evaluation of the sum rules for 1 — 2 derived in 
Sect. 15.21 As the 1 = 2 D-wave irir phase is very small, in 
practice no meaningful constraint results in these partial 



16 



M. Hoferichter, D. R. Phillips, C. Schat: Roy-Steiner equations for 77 — s> ivn 



waves and we therefore restrict ourselves to the S'-wave. 
These sum rules were written explicitly in ((87|) (see also 
App. IC.2|) . The results for and I (2) when the two 
different input tttt phases are chosen are shown in the first 
column of Table [5j The difference between using CCL and 
GKPRY 7T7T phases is very small in both cases. 

Evaluating the sum-rule integrals involves several ap- 
proximations, in particular, we have put Im/io,+ (i) to zero 
above the matching point, neglected partial waves with 
J > 2, and used a resonance approximation for lm fj,±{s). 
Therefore, we must now make sure that the dependence on 
the high-energy part of the integrals and I^ 2 \ higher 
partial waves, and details of the resonance description of 
the left-hand cut, is sufficiently small for the sum-rule con- 
straint to be meaningful. We can estimate the accuracy of 
these approximations by sending the hyperbola parameter 
a — > 00, because in the original setup the dispersion rela- 
tion is independent of a. Thus, any residual dependence on 
a provides a measure of the impact of the approximations 
made. The results of doing this are shown in the second 
column of Table [5j the once-subtracted integral depends 
strongly on the value of a, but the twice-subtracted ver- 
sion is already rather stable under a — > 00. Doubling the 
effect of taking a — > 00 to get a conservative estimate of 
the uncertainty in I^ 2 \ we conclude that 

I (2) = 3.45 ±0.30. (112) 

In order to further test the sensitivity of the sum rule 
to the modeling of the left-hand cut by a set of resonances 
we can check what happens if we switch off resonance 
contributions completely. In the case of 1^ even making 
this crude approximation entails a relatively modest shift 
in the result. As shown in the third column of Table [5J 
the resonances contribute less than 50 % to the full result, 
such that their contribution would have to be drastically 
wrong to exceed the error estimate given in (II 12|) . 

The stability of 1^ under these changes in high-energy 
input makes it worth taking (|112D seriously as a con- 
straint on a particular linear combination of (ai — /?i) /=2 
and (012 — /?2) /=2 - Thus, we will now consider the re- 
sulting sum rule, which arises from the twice-subtracted 
Muskhelishvili-Omncs representation, in more detail. 

First of all, we test if the parameter sets of Table 2] 
fulfill the sum rule. The error analysis is complicated by 
the fact that in the GMM set no uncertainty estimate is 
given for (ai - (5iY , while in the ChPT set the error 
induced by the LEC dependence of (0:2 — P2Y is dif- 
ficult to assess. To obtain a rough estimate, we use the 
ChPT error for («! - fa)* also for GMM, and vice versa 
for (a2 — P2Y ■ This, together with the number (|112p . 
leads to the results summarized in Table [5J We conclude 
that the sum rule is fulfilled for both sets, although rather 
marginally in the case of GMM. (This is mainly due to 
the fact that {012 — P2Y differs quite substantially be- 
tween ChPT and GMM.) The largest uncertainty in the 
sum rule is driven by lack of knowledge of the quadrupole 
polarizability. 



Observing that both dipole polarizabilities as well as 
(c*2 — P2Y have an accurate ChPT prediction, we can 
turn around the argument and use the sum rule to derive 
an improved value for (02 — P2Y ■ Using (|112p and the 
ChPT prediction for the isospin-two dipole polarizability, 
((87)) leads to 

(a 2 - /3 2 ) /=2 = (-18.2 ± 1.3) • lO^fm 5 . (113) 

Using, in addition, the ChPT prediction for (ot2 — foY" > 
we find 

(a 2 - /3 2 ) 7r± = (15.3 ± 3.7) • 10~ 4 fm 5 , 
(a 2 ~ /3 2 ) /= ° = (39.4 ± 6.0) • 10~ 4 fm 5 , (114) 

where the increase in uncertainty compared to (|113[) is due 
to the ChPT uncertainty in (a.2 — P2Y ■ I n the remain- 
der of the paper we will make use of the improved value 
(|114[) when referring to the ChPT predictions for pion 
polarizabilities. Note that, as expected given the results 
of Table [5J our sum- rule value of («2 ~ P2Y is consis- 
tent with the first ChPT number quoted in Table [H but it 
is not consistent with the larger number found when the 
LECs of [63] are taken as input. 

8.2 Total cross section 

Before performing the analytic continuation to the a pole, 
we wish to make sure that the amplitude on the real 
axis is reasonably well described — at least up to y/t — 
1 GeV, which we assess to be the region which will influ- 
ence the analytic continuation to the a pole. The results 
for the cross section are depicted in Figs. [5] and [5J Below 
the matching point, the results for the once- and twice- 
subtracted formulation are provided for both ChPT and 
GMM polarizabilities. The uncertainty due to tttt input, 
represented by the grey band, is estimated by the varia- 
tion between CCL and GKPRY phases and proves to be 
very small. The low-energy region is totally dominated 
by the Born terms in the charged process, but it is very 
sensitive to the a in the neutral reaction. The prediction 
of the twice-subtracted dispersion relation is in especially 
good agreement with 77 — > tt°tt° data (see Fig. [5J, with 
the level of agreement comparable to that obtained in the 
coupled-channel fit of [50] . 

Above the matching point, we exploit the fact that 
the cross section is dominated by the ^(1270), and thus 
can be well approximated by employing a Brcit-Wigncr 
description of this resonance in h,2 = -(t) and putting all 
other partial waves to zero. In this way, (|108[) alone yields 
a good description of the neutral cross section above the 
matching point. In contrast, in the charged case an addi- 
tional background is necessary. As observed in [21], this 
can be most easily achieved by adding the Born terms and 
the off-shell contributions dropped in the transition from 
(I107P to (|108l) back into the charged- channel amplitude for 
t l 2,-(t)- Moreover, after the transition to the isospin basis, 
we add a constant background phase to ensure matching 
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0.4 0.6 0.8 1 1.2 1.4 

s/t I GeV 



Fig. 5. Total cross section for 77 — > tt°tv° 5 , 10 and 77 — > 
tt+tt" [5KUE] for I cos#| < 0.8 and | cos6»| < 0.6, respectively. 




%.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Vt I GeV 



Fig. 6. Total cross section for 77 — > tt°tt° for | cos#| < 0.8 in 
the low-energy region. 




Vt I GeV 



Fig. 7. Dependence of |/i2,_ (t)| on the number of subtrac- 
tions. The grey bands indicate the difference between CCL 
and GKPRY phases. 

with the 7T7T phase below the matching point. However, 
if C? Cl is chosen to be negative, the mismatch of the 
phases is very small: we find a correction of <5 C orr = —0.09 
and S COII = —0.04 in order to obtain agreement with the 
CCL and GKPRY phases, respectively. 

Finally, we comment on the analyticity properties of 
the partial waves at the matching point. As shown in 
the appendix of |48j . the solutions in terms of Omnes 
functions automatically fulfill continuity at the matching 
point, but the derivative at t m is not determined. There- 
fore, in general, strong cusps can occur at the matching 
point. For example, if the background in the charged reac- 
tion is dropped, the neutral cross section above t m is still 
correctly reproduced, but the input for the 1 = compo- 
nent changes, which affects the neutral cross section below 
i m : the result for |/i2^°(0l exhibits a strong cusp below t m , 
which translates into an (unphysical) sharp peak of about 
15 nb in the neutral cross section directly below t m . The 
fact that this effect is much smaller in the full solution pro- 
vides evidence that our model for the high-energy region 
is reasonably accurate, because only a specific input of 7T7t 
phases, polarizabilities, and imaginary parts above t m will 
yield a smooth solution for around t = t m . In the 

language of [SB] such a solution corresponds to an "ana- 
lytic input" . If the input above the matching point were 
sufficiently well known, one could thus derive constraints 
on the polarizabilities by requiring a no-cusp condition. 
These constraints would be similar to those derived in 
[451148] for 7T7t and irK scattering lengths. However, the 
input above the matching point is not very well known in 
77 —¥ 7T7T, so we content ourself with requiring that the 
cusp at t m is not too large, such that the input we are 
using is reasonably close to being "analytic" . 



8.3 Two-photon coupling of the er 



We now present our results for the two-photon width r , (T77 
as a function of the pion polarizabilities. (07 + fix) and 
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4 6 

-A) /=0 / 10- 4 fm 3 

Fig. 8. F a - n as a function of the 1 = pion polarizabilities. 
The black line refers to the unsubtracted case and the colored 
lines to the twice-subtracted version with (a.2 — @2) I=0 as indi- 
cated (in units of 10" 4 fm 5 ). The grey band for the uncertainty 
in the nir input is estimated by the variation found when CCL 
and GKPRY phases and a parameters are chosen. The cross 
corresponds to the twice-subtracted case plus ChPT input. 




4 6 

(ai - A) /=0 / Krfm 3 



10 



Fig. 9. Individual contributions to r ajl for CCL phases and 
a parameters. The dashed lines refer to the once-subtracted 
version. The solid lines refer to the twice-subtracted version, 
with an input value of («2 — P2) I=0 from H114[) . The black lines 
denote the results for a left-hand cut modeled solely by the 
Born terms, while the red, blue, and green lines are obtained 
by successively adding resonances in the limit a — > 00, terms 
for finite a, and D-wave contributions. 



(«2 + h) only feature as subtraction constants in the D- 
waves, which, in turn, influence / n CT77 only indirectly via 
the corresponding coupling to the 5- wave in (I75|) . More- 
over, the imaginary part of these D- waves is dominated 
to a large extent by the /2(1270), and the dependence 
of |/i2^°WI on the number of subtractions is very weak. 
Fig. [3 shows the results when we choose ChPT polari- 
zabilities. The variation between the |^2-°WI solutions 
for different numbers of subtractions is so small that the 
uncertainty in the 7T7t phases (estimated as the difference 





1 subtraction 


2 subtractions 


ChPT 


1.3 ±0.1 


1.7 ±0.4 


GMM 


1.4 ± 0.1 


2.0 ±0.2 



polarizabilities for CCL phases and a parameters. 



between CCL and GKPRY) becomes of comparable size. 
Consequently, the dependence of on (pi\ + /3i) /=0 

and («2 + ^2) I=0 is totally negligible, such that we end up 
with the dipole polarizability (a± — /3i) /=0 as the only free 
parameter that affects r <J ~ fl in the once-subtracted disper- 
sion relation. A second subtraction additionally requires 
the quadrupole polarizability (a2 — f32) I=a as input. The 
resulting correlation between F ail and the pion polari- 
zabilities is depicted in Fig. [3] This model-independent 
correlation is the main result of our study. 

The role of the different contributions to the left- 
hand cut is illustrated in Fig. [9j starting from the Born- 
term approximation (black line) , we add resonances in the 
limit a — > 00 (red line), the additional terms for finite 
a = — 7.5 M% (blue line), and D-wave contributions (green 
line). The twice-subtracted version (solid lines) is hardly 
affected by any of these changes. In the once-subtracted 
case (dashed lines) we see that .D-wave and resonance 
contributions are of comparable size. We therefore expect 
that, as soon as resonances yield a significant contribution 
to the left-hand cut, the coupling between S- and D-waves 
should also become numerically important in any descrip- 
tion of data that is based on a Muskhelishvili-Omnes rep- 
resentation. 

The results for r , (T77 given ChPT and GMM choices for 
the polarizabilities are summarized in Table [71 where the 
errors only include the uncertainties from the pion polari- 
zabilities. We have checked that these numbers are insen- 
sitive to the details of the input above the matching point. 
For the ChPT parameters the results from the once- and 
twice-subtracted equations are consistent. However, there 
is significant tension between these two results in the case 
of GMM. This issue seems to be related to the relatively 
small value of {012 — foY in that polarizability set: in- 
creasing this polarizability and thus bringing it closer to 
the ChPT prediction would both improve the fulfillment 
of the 1 — 2 sum rule (cf. Sect. I8.1|l and bring r aiJ down 
from the (2.0 ± 0.2) keV result given in Table H 

In view of Fig. |H1 we consider the outcome of the twice- 
subtracted system of equations as the most reliable one, 
despite the fact that input for the quadrupole polarizabil- 
ity is required. Although the result based on the GMM 
parameters is ostensibly more precise, the Muskhelishvili- 
Omnes representation used in [30] is, for the reasons ex- 
plained in Sect. 15.41 not fully consistent with the Roy- 
Steiner equations derived here. For this reason, we follow 
the philosophy of 5!) and combine the strict predictions 
of our Roy-Steiner equations with ChPT input for pion 
polarizabilities to obtain our final result 



(1.7±0.4)keV, 



(115) 
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Reference 






Penn. 06 




4.1 ±0.3 


MNO 08 


m 


3.9 ±0.6 


T\ TTVTTX T 1 1 

JVUN W 1 1 


U2J 


3.08 ± 0.82 


FK 06 


US] 


0.62 


Penn. et al. 08 


[26]. sol B 


2.4 ±0.4 


ORS 08 


m 


1.8 ±0.4 


OR 08 


m 


1.68 ±0.15 


BP 08 


\m 


1.20 ±0.40 


Mao et al. 09 


m\ 


2.08 



Table 8. Previous results for r , (T77 in keV. 



Penn. 06 - I 1 1 

MNO 08 - I 1 1 

MNW11- I 1 1 

FK 06 - + 

Penn. et al. 08 - I 1 1 

ORS 08 - I 1 1 

OR 08 - I — I — I 

BP 08 - I 1 1 

Mao et al. 09 - + 
This work - I 1 1 

1 2 3 4 5 

IV,. / keV 

Fig. 10. Previous results for r all with acronyms defined in 
Table [8] The red point corresponds to our result. 

which is depicted by the cross in Fig. [5] We note that 
the uncertainty here is broad enough to encompass all 
the central values in Table [7] A comparison with previous 
results for r , (T77 is shown in Tableland Fig. [TU1 

9 Summary and conclusion 

In this paper we constructed a full set of Roy-Steiner 
equations for 771 — > 771. In particular, we worked out all 
necessary integral kernels for zero, one, and two subtrac- 
tions explicitly up to D-waves, and directly identified the 
subtraction constants with the pion polarizabilities. We 
studied the range of validity of the system, and found 
that the equations for the crossed channel 77 — > mr are 
rigorously valid up to 1 GeV — a domain that comfortably 
includes the a pole. Extending this range requires addi- 
tional assumptions. Truncating the system at J = 2, we 
then concentrated on the equations for 77 — > 7r7r, whose 
solution in terms of a Muskhelishvili-Omnes representa- 
tion with a finite matching point was discussed. Compar- 
ing our equations with existing approaches in the litera- 
ture, we found a coupling between S- and D-waves, which 



has previously been neglected in calculations of 77 — > irir 
based on dispersion relations, but seems to be numerically 
comparable to the contributions of resonances in 77 — > irn 
to the left-hand cut. 

Demanding continuity of the Muskhelishvili-Omnes 
representation at the matching point, we derived sum rules 
for the 1 = 2 partial waves that relate dipole and quadru- 
pole polarizabilities to integrals over the left-hand cut. We 
used the S'-wave sum rule, together with ChPT input for 
the neutral-pion quadrupole polarizability and the dipole 
polarizabilities, to obtain a new, more accurate, prediction 
for the charged-pion quadrupole polarizability 

(a 2 ~ /3 2 ) w± = (15.3 ± 3.7) • l(T 4 fm 5 . (116) 

While the central value is hardly shifted, the error esti- 
mate is difficult to obtain in ChPT alone due to a strong 
dependence on poorly known low-energy constants. 

The main application of our formalism though, was 
a study of the two-photon width of the a resonance. To 
this end, we first showed that the cross section for both the 
charged and neutral channel can be accurately reproduced 
by approximating the high-energy region above 1 GeV by 
a Breit-Wigner ansatz for the /2(1270) resonance and em- 
ploying a suitably chosen background amplitude in the 
charged case. With this input, our Muskhelishvili-Omnes 
representation yields a good description of the available 
data in the low-energy region. 

We then presented the results of the analytic continu- 
ation to the a pole as a correlation plot between the per- 
tinent pion polarizabilities and r a ~ n . Our most general 
finding is a correlation between 7 = pion polarizabilities 
and r , - 77 . This correlation results solely from analyticity, 
unitarity, crossing symmetry, and the accurately known 
7T7r phases in the region below 1 GeV. We also provided 
a specific result for r ai7 by taking ChPT predictions for 
the pertinent pion polarizabilities. We therefore look for- 
ward to the results of the ongoing efforts at COMPASS to 
extract pion polarizabilities, which — as shown here — will 
further improve knowledge of r ail . 
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A Kernel functions for the s-channel projection 

The functions Nf(s) for J < 2 are given by 

s M 2 r 1 -1 

NUs) = MfS ~ "^K ^o(y) - 2} + AN t^ 

M 2 M 2 r 1 -1 

1 M 2 M 2 r 1 i 

N * {8) = ~3Mf~ s + -r{ d l-i(y)Qo(y) g + y) + an 2 -(s), 

s + M 2 



where 



and 



Qo(z) = ~ 



da; 

2 J z-x'' 
1 



s-MV 



Qo(z±ie) = - log 



l + z 



1 - z 



Ti^{l-z% 



denotes the lowest Legendre polynomial of the second kind. The remainders ANf (s) 



AN+ {0 \s) = 0, AN? (1 >(s) 



3M n a 



(ai + A)q 2 , ^^(s) = Z\A^(s) - 



lSikf^a 



(a 2 + /3 2 )q 4 



Z\iV 2 +(0) (s) = Z\W 2 + w (s) = 0, AN^ A '(s) = 



+ ( 2 )^ - 



SOikf^a 



(a 2 + /3 2 )q 4 



zi7V- (0) (s) 
Z\A^ (2) (s) 



2q 2 



3(M*-a) : 



1 



3a 



(ai+ftjq 4 , 



2M T 
3a 



(ai-/3!)q 2 -^(a 2 - (,)q' 



(a 2 +/3 2 )q 6 , 



z\7v 2 - (0) ( s )-o, z\tv 2 - (1) ( s ) = 



Toh^ + ^ (2) ^ = wa^ ~ ^ q4 + idb 



(A.1) 
(A.2) 

(A.3) 



(A.4) 



(a 2 + /3 2 )q 6 



contain the pion polarizabilities according to the number of subtractions indicated by the superscript. 
A.l s-channel 

The kernels for the unsubtracted case read 



K++(s,s') 
K+ + (s,s') 



sq I 1 1 3 

s'q' 2 1 s' — s s' — a 2q 2 

* 2 ( ' 1 1 \ 5 



: Qo(x s ) — 



■sq 

s'q' 2 IV s' — s s' — a J 3 



(2/3 + a-l) 



5 
2q^ 



(2a£ - l)Q (z a ) - (J« + /3- 5 



a+2 -Jj^- 2 x: 



s 3/ (l + x s ) 2 /o 1 3z s a; 2 

'^2l^^ ( s " )0o(a;s) " 6 " T - T 



1 



1 



s q^ 

s'q' 2 1 \ .-. s • 



2q 2 



(l + x s ) 2 



(2z 5 - 1)(2^ - l)Q (x s ) 



-(!-5)-(f + i*-i)«.-(^|-5)*5— ! 
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'> = 2^{ (7b - 7^) (^ 2 - 2^ - 2 ^ + 2 - 2^) }■ 

K * +{a > s ' } = { {-^h vh) { ( 5/3 - 3a - D q2 + ir ( 8a2 + 5a(3 - Afj) + 10 ^ - 3/3 + 2/ ^ I 



+ 



5 s' - s 



X 2 -+( S , S ') 



K^ + (a,a') 



6 s' — a 

2 / ! 



2s'q' 2 Vs'-s s'-a/lO 



(q 2 -q' 2 «), 



-^{(Fb-F^)((^T-^-<C^-^)-fF^}- 



*r 2 - (*>*') 



1 



q' 2 I s' — s s' — a 2q 2 



-. Qo{x s ) + — — 



1 



(l-^ s ) 2 



(24 + l)Qo(^) 



1 



a -+- - \ x s - —x. 



1 3 J (l-z s ) 2 . . 1 3x s x 2 s 1 



1 



1 



(1 - x s 



q' 2 I \ s' — s s' — aj 2q 2 

f/3 1\ fa 3 3\ / 3 1\ 2 
(3 + 6)-(3-2 /3 -l)^-( /3 -2 a+ 2)^- 



(2x s + l)(2x' s + l)Q (x s ) 



ax" 



(A.5) 



where 



a = 



q" s — a 
q' 2 s' — a ' 



/3 = l-a- S ' " S + S 1 /2 2M2 ; 4 = «^ + A ^ = 1- S + S L 2 2M " . < = "^ + /3. (A.6) 



s'-a 2q' 2 ' s s ^' s 2q 2 
The corresponding versions for the once-^ 1 )) and twice-(( 2; ) subtracted case can be obtained by adding 

__(2J + !)-*„ I- 

2s q 4 



AK+p 1 \s,s') = -^(2J' + 1)^- 
AK^\s,s') = ^(2J' + l)(^- 



5 



s' -Ml s'-aj l + z' a 
2 



1 \g t d J n(z' s ) 



3' -Ml s'-aj 1 + z 
1 \d{ 1 (z' s ) 



1 d{[{z' s ) 



(s' - Af 2 ) 2 1 + 
/ 1 n2 



^7;( 2 ) (s , s ') = -^(2j' + i)< 



q 



15 J 105 J 



1 ^/u«) 



s' — M 2 s' - a J 1 + ^ 
2 

s'-M; 
2 



) -4M 2 (V 

1 dl'M) 



ji 



(«' - M 2 ) 2 1 + 



_ _ _i_) a<( , _ 4M3 fg|| o - <^W«' - ^iSfl J}. 



(2) 2q 4 / ^2\// 2 1 ^ a 



(,s'-M 2 ) 2 l-z' s 



(A.7) 
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with 



,, s' + a ,. _ Ml -a 

I n , 1 Gt^s In 



s\0 



s' — a 



" s| ° 2q' 2 (s'-a)' 



A. 2 t-channel 



The non- vanishing unsubtracted kernel functions are 

G^T (*>*') = 2t , ( ^4 M 2 ) {( 1 +^) 2 ( 2 ^ - VQoixt) - I - 3xt - 2xl), 

Gl + (s,t') = -^{(1 - x*) 2 Q (z t ) + 2 - st}, 

G 2 - + (s,t') = -^{(1 - x t ) 2 (2x t + l)Q (x t ) -^+3x t - 2x 2 t y 

G- l2 + (s,t') = ^{(1- ^) 2 P2(x;)Qo(x t ) - (1 + 2 7 - 36) + (± + 3 7 - ™)x t - ^f}, 
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5V6q 2 
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1 ^22 
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2q 
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The contributions from Born terms and subtraction constants for J < 2 are 



2<r(i) J ' 



B.l s-channel 

The non-vanishing kernel functions for the unsubtracted case are 

3i r , 2 37 + 3(5- 1 
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The non- vanishing kernel functions for J, J' < 2 arc 
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C Omnes solutions for the 77 — >• 7T7t partial waves 
C.l / = 

The equations hold for < <5(t m ) < w, which all the 1 = partial waves satisfy. We only show the results for the once- 
and twice-subtracted version, as the reduction to the unsubtracted case is straightforward 

hoAt)=^ + (t) + —(a 1 -p 1 )w (t) + ——i j dt ^^^ +ydt _ t)mn , 

h . + (t) = igj.(t) + - ft)t(l - * rt (0))/?o(t) + ^("2 - I32)t 2 f2 (t) 



M. Hoferichter, D. R. Phillips, C. Schat: Roy-Steiner equations for 77 — s> ivn 



25 



+ 



t 3 n (t) 



v; I J t«(f -t)|/? (*')| J t' 3 {t'-t)\Q Q {t')\ 



ft2,+ (t) - 4 2 ,+ (t) + 



t 2 (t-4Af2)!2 2 (t) 



dt'- 



sin5 2 (i')^2,+(t') 



i'2( t '_4M|)(t'-t)|/2 2 (i')| 7 i' 2 (t'-4Af2)(t'-<)|f2 2 (t')| 



dt'- 



lmh 2 ,+{t') 



h 2 ,.(t) = A£Ut) + * ( * J^\ a! + P l )f2 2 {t) 



10y/GM w a 



t 2 (t-4M*)f2 2 (t) 



dt' 



sm8 2 {t')A { 2 l) _{t') 



t> 2 (t>-4M%)(t>-t)\f2 2 {t>)\ J t' 2 (t> -4M%)(t> -t)\f2 2 (t')\ 



At' 



lmh 2 ^(t') 



h 2 ,.(t) = A^_(t) + 



w 2) tft - 4M 2 ) 



10\/6M w a 



(ai + - t fi 2 (0)) + — (a 2 + ft)] « 2 (t) 



t 3 (t-4M2)/? 2 (t) 



dt'- 



i' 3 (t'-4M2)(t'-t)|tt 2 (i')| 7 t'3(t> -4M%)(t> -t)\f2 2 (t>)\ 

K 4Ml tm 



dt'" 



Im/i 2> _(t') 



(C.l) 



The left-hand-cut functions Aj^ ± (t) are defined in (|75j) and 



C.2 J = 2 

The solutions for — 7r < 5(t m ) < 0, which is true for all I = 2 partial waves, read 



V+(*) = <|(i) + ^(ai-A)^o(t) ( 



en (t)(t m - 1) j 7 dt , smSo(t')A^(t') 



MM2 t m ' 

/*,,+ (*) = A^(t) + ^{a x - fr)t\l + t m A,(0))1 rtoW^^ + ^(«2 - fe)^^^ 



. Tmho, + (t') 

\j t^(t m -t')(t' -t)\Q (t')\ ' J t' 3 (t m -t')(t' -t)\Q (t')\ 

AMI *m 



h2 . +{ t ) = A 2 , + it ) + t ^^mm^A\ , (l , 



smd 2 {t')A 2<+ (t') 



4i\/H 



t' 2 (t'-4MZ)(t m -t>)(t>-t)\n 2 (t')\ 



dt' 



Im h 2 , + (t') 



ti*(t> - 4M*){t m -t>)(t> -t)\n 2 {t>)\ 



h 2 ,_(t) = AP_(t) + ^ J^V + ftW)-*" 1 ' 



t 2 {t-AMl)[2 2 {t){t m -t) 



dt' 



sin5 2 (t')^ 2 1 -(*') 



dt' 



lmh 2 ,-{t') 



t> 2 (t> -AM 2 ){t m -t')(t' -t)\n 2 {t')\ 



t> 2 (t> -4M*){t m -t')(t' -t)\Q 2 (t')\ 



h 2 .At) = ^_(t) + t -^-0^.[{ ai + h) [1 + £(1 - t m n 2 (o))] + ^(a 2 + p 2 )}n 2 {t) 



Cm t 
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t 3 {t-4M^)Q 2 (t)(t m -t) 



dt' 



wi8 2 (t')Af_{t') 



MAP 



t'Ht* -AM%){t m -t')(t' -t)\n 2 {t')\ 



dt' 



Imh 2 ,-( t ') 



t' 3 (t'-4M^)(t m -t')(t'-t)\Q 2 (t')\ 



(C.2) 



The sum rules discussed in Sect. 15.21 are then 



= 



= 





(ai - 


2a 


AU 


(ai - 




2a 


t 3 


' tm 

/ 


IT 1 





tr 



7T 



, sinJoCtO^g+C*') , 7 , + , Imho,+(*') 



dt 



dt' 



t' 2 (t'-t m )|a,(t')| ' J t>z(t> -t m )\f2 (t>)\ j' 

M 

^o(0)) + ^(a 2 -^ 2 )4 

t(2 



, sinM*'K%(i') , /-^ Imft , + (t') 



f^t'-t^W)! J t' 3 (t'-t m )\n (t')\ 



f^(t m -4M 2 ) 



dt'- 



smS 2 {t')A 2<+ (t') 



f' 2 (t'-4M 2 )(t'-t m )W)l 
MM2 t 



= p= [ai + Pi) 



4(t m -4M2) 



t m (t m ~ 4M 2 ) 

ioV6Af„.a 

4(t m -4M 2 ) 



dt'- 



t' 2 (t'-4M 2 )(t'-i m )|r? 2 (t')| 

"4 A/2 « 



(ai + - t m fl 2 (0)) + y|(« 2 + A) 



dt'- 



sin^t')^ 2 ^') 



dt' 



Im/i 2 ,+ (t') 



t' 2 (t'-4M 2 )(t'-t m )|tf 2 (t')| [' 



dt' 



Im/i 2 ,_(t') 



t' 2 (t'-4M2)(t'-t m )|/2 2 (t')| (' 



dt'- 



Im h 2 -(f) 



4 a/2 



f3(f-4Af2)(f -t m )|fl 2 (f)| '7 t' 3 (t'-4Af 2 )(t'-t m )|^ 2 (t')| 



(C.3) 



References 

1. I. Caprini, G. Colangelo and H. Leutwyler, Phys. Rev. 
Lett. 96 (2006) 132001 [arXiv:hep-ph/0512364| . 

2. S. M. Roy, Phys. Lett. B 36 (1971) 353. 

3. S. Weinberg, Physica A 96 (1979) 327. 

4. J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142, 
Nucl. Phys. B 250 (1985) 465. 

5. H. Marsiske et al. [Crystal Ball Collaboration], Phys. Rev. 
D 41 (1990) 3324. 

6. J. Boyer et al., Phys. Rev. D 42 (1990) 1350. 

7. H. J. Behrend et al. [CELLO Collaboration], Z. Phys. C 
56 (1992) 381. 

8. T. Mori et al. [Belle Collaboration], J. Phys. Soc. Jap. 76 

(2007) 074102 |arXiv:0704.3"538l [hep-ex]]. 

9. S. Uehara et al. [Belle Collaboration], Phys. Rev. D 78 

(2008) 052004 arXiv:0805.3387 [hep-ex]]. 

10. S. Uehara et al. [BELLE Collaboration], Phys. Rev. D 79 

(2009) 052009 |arXiv:0903.3697l [hep-ex]]. 

11. A. I. L'vov, V. A. Petrun'kin and M. Schumacher, Phys. 
Rev. C 55 (1997) 359. 

12. D. Drechsel, B. Pasquini and M. Vanderhaeghen, Phys. 
Rept. 378 (2003) 99 [arXiv:hep-ph/0212124| . 



13. J. Bernabeu and J. Prades, Phys. Rev. Lett. 100 (2008) 
241804 :i arXiv:0802.1830l [hep-ph]]. 

14. M. Schumacher, Eur. Phys. J. C 67 (2010) 283 
|arXiv:1001.0500l [hep-ph]] . 

15. M. R. Pennington, Phys. Rev. Lett. 97 (2006) 011601. 

16. G. Mennessier, S. Narison and W. Ochs, Phys. Lett. B 665 
(2008) 205 |arXiv:0804.4452l [hep-ph]]. 

17. G. Mennessier, S. Narison and X. G. Wang, Phys. Lett. B 
696 (2011) 40 [arXiv: 1009.27731 [hep-ph]]. 

18. N. Muskhelishvili, Singular Integral Equations, P. Noord- 
hof, Groningen, 1953. 

19. R. Omnes, Nuovo Cim. 8 (1958) 316. 

20. M. Gourdin, A. Martin, Nuovo Cimento 17 (1960) 224. 

21. O. Babelon, J. L. Basdevant, D. Caillerie, M. Gourdin and 
G. Mennessier, Nucl. Phys. B 114 (1976) 252. 

22. D. Morgan and M. R. Pennington, Z. Phys. C 37 (1988) 
431 [Erratum-ibid. C 39 (1988) 590], Phys. Lett. B 272 
(1991) 134. 

23. J. F. Donoghue a nd B. R. Ho lstein, Phys. Rev. D 48 (1993) 
137 [arXiv"hep- ph/9302203 . 

24. D. Drechsel, M. Gorchtein, B. Pasquini and M. Van- 
derhaeghen, Phys. Rev. C 61 (1999) 015204 
[arXiv:hep-ph/9904290] . 



M. Hoferichter, D. R. Phillips, C. Schat: Roy-Steiner equations for 77 — s> ivn 



27 



25. L. V. Fil'kov and V. L . Kaslievarov, Phys. Rev. C 73 
(2006) 035210 [arXiv:nucl-th/0512047] . 

26. M. R. Pennington, T. Mori, S. Uehara and Y. Watanabe, 
Eur. Phys. J. C 56 (2008) 1 |arXiv:0803.3389l [hep-ph]]. 

27. J. A . Oiler, L. Roca an d C. Schat, Phys. Lett. B 659 (2008) 
201 larXiv:0708.f659l [hep-ph]]. 

28. J. A. Oiler and L. Roca, Eur. Phys. J. A 37 (2008) 15 
|arXiv:0804.0309l [hep-ph]] . 

29. Y. Mao, X. G. Wang, O. Zhang, H. Q. Zheng 
and Z. Y. Zho u, Phys. Rev. D 79 (2009) 116008 
|arXiv:0904.f445l [hep-ph]] . 

30. R. Garcfa-Martm and B. Moussallam, Eur. Phys. J. C 70 
(2010) 155 [arXiv:1006.5"373l [hep-ph]]. 

31. Yu. M. Antipov et al, Phys. Lett. B 121 (1983) 445. 

32. T.A. Aybergenov et al, Sov. Phys.-Lebedev Inst. Rep. 6 
(1984) 32; Czech. J. Phys. B 36 (1986) 948. 

33. J. Ahrens et al., Eu r. Phys. J. A 23 (2005) 113 
[arXiv:nucl-ex/04070TT] . 

34. J. Bijnens and F. Cornet, Nucl. Phys. B 296 (1988) 557. 

35. J. F. Donoghue, B. R. Holstein and Y. C. Lin, Phys. Rev. 
D 37 (1988) 2423. 

36. S. Bellucci, J. Gasser and M. E. Sainio, Nucl. Phys. 
B 423 (199 4) 80 [E rratum-ibid. B 431 (1994) 413] 
[arXiv:hep-ph/9401206j . 

37. U. Btirgi, Phys. Lett. B 377 (1996) 147 
[arXiv:he p-ph/9602421 ; Nucl. Phys. B 479 (1996) 
392 arXiv:hep-ph/9602429 . 

38. J. Gasser, M. A. Ivanov and M. E. Sa inio, Nucl. Phys. B 
728 (2005) 31 [arXiv:hep-p h /0506265] . 

39. J. Gasser, M. A. Ivanov and M. E. Sa inio, Nucl. Phys. B 
745 (2006) 84 [arXiv:hep-ph/0602234] . 

40. A. V. Guskov, Phys. Part. Nucl. Lett. 7 (2010) 192. 

41. A. Guskov [COMPASS Collaboration], Nucl. Phys. Proc. 
Suppl. 198 (2010) 112. 

42. A. R. Edmonds, "Angular momentum in quantum me- 
chanics", Princeton University Press, Princeton (1960). 

43. M. Jacob and G. C. Wick, Annals Phys. 7 (1959) 404 [An- 
nals Phys. 281 (2000) 774]. 



44. I. Guiasu and E. E. Radescu, Annals Phys. 120 (1979) 
145. 

45. B. Ananthanarayan, G. Colangelo, J. Gasser 
and H. Leutwyler, Phys. Rept. 353 (2001) 207 
[arXiv:hep -ph/0005297] . 

46. G. E. Hite and F. Steiner, Nuovo Cim. A 18 (1973) 237. 

47. B. Ananthanarayan and P. Biittiker, Eur. Phys. J. C 19 
(2001) 517 [arXiv:hep-ph/0012023l . 

48. P. Biittiker, S. Descotes-Genon and B. Moussallam, Eur. 
Phys. J. C 33 (2004) 409 [arXiv:hep-ph/0310283j . 

49. F. E. Low, Phys. Rev. 96 (1954) 1428. 

50. S. Mandelstam, Phys. Rev. 112 (1958) 1344. 

51. G. Hohler, in Landolt-Bornstein, 9b2, ed. H. Schopper, 
Springer, Berlin, 1983. 

52. C. Itzykson and J.-B. Zuber, Quantum Field Theory, Mc 
Graw-Hill Inc, New York (1980). 

53. H. Lehmann, Nuovo Cim. 10 (1958) 579. 

54. K. M. Watson, Phys. Rev. 95 (1954) 228. 

55. L. Epele and G. Wanders, Phys. Lett. B 72 (1978) 390; 
Nucl. Phys. B 137 (1978) 521. 

56. J . Gasser and G. Wand ers, Eur. Phys. J. C 10 (1999) 159 
[arXiv:hep-ph/9903443l . 

57. G. Wanders, Eur. Phys. J. C 17 (2000) 323 
[arXiv:hep -ph/0005042] . 

58. I. Caprini, G. Colangelo, H. Leutwyler, in preparation. 

59. G. Colangelo, J. Gasser and H. Leutwyler, Nucl. Phys. B 
603 (2001) 125 [arXiv:hep-ph/0103088j . 

60. R. Garcfa-Martm, R. Kaminski, J. R. Pelaez, J. Ruiz de 
Elvira and F. J. Y ndurain, Phys. Rev. D 83 (2011) 074004 
|arXiv:1102. 21831 [hep-ph]]. 

61. I. Caprini, private communication. 

62. R. Garcfa-Martm, R. Kaminski, J. R. Pelaez and 
J. Ruiz de Elvira, Phys. Rev. Lett. 107 (2011) 072001 
|arXiv: 1107. 16351 [hep-ph]] . 

63. J . Bijnens and J. Prad es, Nucl. Phys. B 490 (1997) 239 
[arXiv:hep-ph/9610360l . 

64. K. Nakamura et al. [Particle Data Group], J. Phys. G 37 
(2010) 075021. 



